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Landsat-4 Scientific Characterization 
Early Results Symposium 

Tuesday, 22 February 

Bldg. 8, Auditorium 

7:30 a.m. Registration 

9:30 a.m. First General Session 

Introduction, Description and Highlights 

12:00 noon Lunch 

1:30 p.m. Second General Session 

Representative Early Results 


5:30 p.m. 


Adjourn 


Landsat-4 Scientific Characterization 
Early Results Symposium 

Wednesday, 23 February 

Bldg. 8, Auditorium 

7:30 a.m. Registration 

8:00 a.m. Third General Session 

Scientific Characterization 

9:30 a.m. Adjourn 

10:00 a.m. TM Introductory Session 

12:30 p.m. Lunch 

1:30 p.m. TM Radiometry and Geometry Session 

5:20 p.m. Adjourn 


Bldg. 3, Auditorium 

10:0^ a.m. MSS Session 

MSS Investigations Results 

12:30 p.m. Lunch 

1:30 p.m. MSS Investigations Results, Continued 

2:30 p.m. MSS Methods and Results 

5:30 p.m. Adjourn 


Landsat-4 Scientific Characterization 
Early Results Symposium 


8:00 a.m. 
8:45 a.m. 
12:00 Noon 
1:10 p.m. 
3:30 p.m. 
4:00 p.m. 

5:00 p.m. 


8:30 a.m. 
12:10 p.m. 
1:30 p.m. 
3:30 p.m. 


Thursday, 24 February 
Bldg. 8, Auditorium 

Registration 

TM Radiometry and Geometry Session, Continued 
Lunch 

TM Radiometry and Geometry Session, Continued 
Adjourn 

Fourth General Session 

Research Plans and Panel Discussion 

Adjourn 


Bldg. 3, Auditorium 

TM Applications Session 
Lunch 

TM Applications Sessions, Continued 
Adjourn 



10:20 a.m. 


Welcome and Introduction 

Vince Salomonson 

Landsai-4 Project Scientist, Goddard Space Flight Center 

Landsal-4 Program 

William Webb 

Landsat-4 Mission Operations Manager, Goddard Space Flight Center 

NOAA'a Landsal Program 

Russell Koffler 

Director, Office of Satellite Data Processing and Distribution, National Oceanic and 
Atmospheric Administration 


10:30 a.m. 


Landsal-4 System Description 


Ted Aepli 

Manager, Landsat Systems Engineering, General Electric 

11:00 a.m. Instrument Descriptions 

Jack Engel 

Manager, Thematic Mapper Program. Santa Barbara Research Center 

11:30 a.m. Landsat-4 Highlights 

Darrel Williams 

Landsai-4 Assistant Project Scientist, Goddard Space Flight Center 
12:00 Noon Adjourn 
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1:30 p.m. 

2:00 p.m. 

2:30 p.m. 

3:00 p.m. 

3:30 p.m. 
4:00 p.m. 

4:30 p.m. 
5:00 p.m. 
5:30 p.m. 


Lftndsat-4 Scientific Characterization 
Early Results Symposium 

Second General Session 
Representative Early Results 
Tuesday, 22 February 

Bldg. I, Auditorium 
John Barker, Moderator 


Earth Imaging and Data Processing for Mapping and Analysts — A Technological 
Assessment 

Ralph Bernstein 
IBM Corporation 

Interpretability of California’s Agricultural Resources on Landsal-4 TM and MSS 
Imagery — An Early Appraisal 

Stephen DeGloria 

University of California, Berkeley 

Evaluation of TM Performance as Applied to Hydrocarbon Exploration 

Jon Dykstra 

Earth Satellite Corporation 

Geologic Utility of l.andsal-4 TM Data — Death Valley, California 

Mike Abrams 

NASA. Jet Propulsion Laboratory 

Break 

An In.tial Analysis of Landsal-4 Thematic Mapper Data for the Discr’ininalion of 
Agricultural Forested Wetland and Urban Land Cover 

Dale Quattrochi 

NASA, Earth Resources Laboratory 

Canrdian Plans for TM Data 

W. Murray Strome 

Canada Centre for Remote Sensing 

Landsa‘-4 Image Data Analysis 

Robert MacDonald 
NASA, Johnson Space Center 

Adjourn 
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8:90 «.m. 


8' IS i.m. 


8:30 i.m. 


9:00 i.m. 


Landsat-4 Scientific Characterization 
Early Results Symposium 

Third General Session 
Scientific Characterization 
Wednesday, 23 February 

Bldg. 8, Auditorium 


landsat-4 Evaluation Program and Scientific Characterization Activities 

John Barker 

Landsat-4 Associate Project Scientist, Goddard Space Flight Center 

MSS Report 

William Alford 

Landsat-4 Senior Systems Programmer, Goddard Space Flight Center 

TM Report 

John Barker 

Landsat-4 Associate Project Scientist. Goddatd Space Flight Center 

Analysis of TM Data from Geologic Terrain in a California- Nevada Scene 

Nicholas Short 

Goddard Space Flight Center 

Rupert Haydn 

Ventralstelle fur Geo-Photogrammetrie w. Ferncrkundung, Federal Republic of 
Germany 


9:30 a.m. 


Adjourn 



Landsat-4 Scientific Characterization 
Early Results Symposium 


10:00 *.m. 


10:30 a.m. 


11:00 a m. 


11:30 a.m. 


12:00 Noon 


TM introductory Session 
Wednesday , 23 February 

Bldg. S, Auditorium 
John Barker, Moderator 


Sensor Radiometry 

J. Barker 

Goddard Space Flight Center 

Sensor Geometry 
J. Engel 

Santa Barbara Research Center 

Radii metric Processing 
J. Barker 

Goddard Space Flight Center 

Geometric Processing 

E. Beyer 
General Electric 

Image Data Quality 

J. Lyon 

Goddard Space Flight Center 


12:30 p.m. 


Adjourn 
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Landsat-4 Scientific Characterization 
Early Results Symposium 

MSS Session 
Wf^mdiy, 23 February 

Bldg. 3, Auditorium 
William Alford. Moderator 

MSS Investigations Results 

MSS Radiometric Accuracy Assessment 

W. Alford 

Goddard Space Right Center 

Spectral Characterization of the Landsat-4 MSS Sciuori 

B.L. Markham 

Goddard Space Right Center 

Investigation of Radiometric Properties of Landial-4 MSS 

W.A. Malila 

Environmental Research Institute of Michigan 

Evr.iuaUoa of Landsat-4 Multispectral Scanner Ground Segment Performance 

A. Zobrist 

. NASA, Jet Propulsion Laboratory 

Radiometric Calibration and Geocod ed Precision Processing of LandsaM MSS 
Prc icts by the Canada Centre for Remote Sensing 

J. Murphy 

Canada Centre for Remote Sensing 

Landsat-4 Multispectral Scanner Radiometric Intraband Performance 

E. Eliason 

U.S. Geological Survey, Ragstaft 

Landsat MSS Scene-to-Scene Registration Assessment 

J.E. Anderson 

NASA, Earth Resources Laboratory 

Geometric Accuracy of Landsat-4 Image Data 
R. Welch 

Department of Geography, University of Georgia 

Geodetic Accuracy of Landsat-4 MSS Data 

L. Oleson 

U.S. Geological Survey, EROS Data Center 



12:30 p. 


Lunch 


Landsat-4 Scientific Characterization 
Early Results Symposium 

MSS Session, Continued 
Wednesday, 23 February 

Bldg. 3, Auditorium 
William Alford, Moderator 

MSS Investigations Results 

1:30 p.m. MSS Geometric Accuracy Assessment 

M. Imhoff 

Goddard Space Right Center 

1:45 p.m. MSS Geometric Performance: Llne-to-LInc Displacement Analysis 

L. Fusco 

Earthnet Program Office, Frascati, Italy 

2:00 p.m. Impact of Landsat MSS Sensor Differences on Change Detection Analysis 

W. Likens 

NASA. Ames Research Center 

2:15 p.m. Comparison of Land Cover Information from Landsrt MSS aud Airborne TMS for 
Hydrological Studies 

J. Oervin 

Goddard Space Flight Center 

MSS Methods and Results 

2:30 p.m. Landsat-4 MSS Radiometric Correction 

A. Singh 
General Electric 

3:30 p.m. Break 

4:00 p.m. Landsat-4 MSS Geometric Correction 

J. Brooks 
General Electric 



5:30 p.m. 


Adjourn 



L, ^dsat-4 Scientific Characterization 
Early Results Symposium 

TM Kudiometry and Geometry Session 
Wednesday, 23 February 

Bldg. I, Auditorium 
Jolia Barker, Moderator 

1:30 p.m. Spectral Characterization of the Landsat I hematic Mapper Season 

B.L. Markham 

Goddard Space Flight Center 

1:90 p.at. Prelaunch Absolute Radiometric Calibratioa of Reflective Beads oa Laadsat-4 

Photcnight Thematic Mapper 

J. Barker 

Goddard Space Flight Center 

2:09 p.m, Characterisation of Radiometric Calibratioa of Laadsat-4 Thematic Mapper Reflective 
Bands 

J. Barker 

Goddard Space Flight Center 

2:25 p.m. Evelaatioa of the Radiometric integrity of Laadsal-4 Band 6 Data 

J. Schott 

Rochester School of Technology 

2:49 p.m. Break 

3:19 p.m. Thermal Band Characterization of Laadsal-4 Thematic Mapper 

J. Lansing 

Santa Barbara Research Center 

3:40 p.m. A Preliminary Assessment of Landaat-4 TM Data Processing by NASA 

D Goodenough 

Canada Centre for Remote Sensing 

4:00 p.m. Preliminary Evaluation of Radiometric Calibration of Laadiat-4 TM Data by the 
Canada Centre for Remote Sensing 

J. Murphy 

Canada Centre for Remote Sensing 

4:20 p.m. Evaluation of the Radiometric Quality of the TM Data Jslng Clustering and 
Multlspectral Distance Measures 

L. Bartolucci 
Purdue University 

4:40 p.m. Analysis and Evaluation of the Landsat-4 MSS and TM Sensors and Ground Data 
V rocessiug Systems — Early Results 

R. Bernstein 
IBM Corporation 


9:20 p.m. 


Adjourn 


Landsat-4 Scientific Characterization 
Early Results Symposium 


1:45 a.m. 

9:00 a. a. 

9:20 a.m. 

•c a.m. 

10:00 a.m. 
10:30 a.m. 

10:50 a.m. 

11:10 a.m. 

11:30 a.m. 


TM Radiometrv sr.u Geometry Session, Cont. 
Thursday, 24 February 

Bldg. I, Auditorium 
John Barker, Moderator 


la Progress Report: Absolute In-Flight C alibration of Landsal-4 Sensors 

P. Slater 

University of Arizona, Optical Science Center 

Landsat-4 TM Calibration and Atmospheric Correction 

W. Hovis 
NOAA/NESDIS 

Scan Angle and Detector Effects in TM Radiometry 

M. Metzler 

Environmental Research Institute of Michigan 

TM Spectral Dimensionality and Data Structure 

E. Crist 

Environmental Research Institute of Michigan 

Break 

MTF Analysis of Landsat-4 TM 

R. Schowengerdt 
University of Arizona 

Thematic Mapper Intraband Radiometric Performance 

H. Kieffer 

U.S. Geological Survey, Flagstaff 

Evaluation of Radiometric and Geometric Characteristics of l.andsat-l TM 

M. Podwysocki 
U.S. Geological Survey 

JSC TM Scrounge Data Preprocessor 

J. Gilbert 

NASA, Johnson Space Center 


I ^ndsat-4 Scientific Characterization 
Early Results Symposium 


11:50 a.m. 

12:00 Noon 
1:10 p.m. 

1:30 p.m. 
2:00 p.m. 

2:20 p.m. 
2:40 p.m. 
2:50 p.m. 


TM Kadiometry and Geometry Session, Cont. 
Thursday, 24 February 

Bldg. 8, Auditorium 
John Barker, Moderator 

Processing of Thematic Mapper Simulator l>ala 

S. Unger 

Goddard Space Flight Center 

Lunch 


r 

H ' 








Evaluation of Lands:.. -4 Thematic Mapper Ground Segment Performance Without 
Ground Control 

A. Zobrist 

NASA, Jc. Propulsion Laboratory 


Early Evaluation of TM Geometric Correction System 

E. Beyer 
General Electric 


The Use of Linear Feature Detection to Investigate TM Data Performance and 
Processing 

C. Gurney 

Systems & Applied Sciences Corp. 

Spatial Resolution Estimation of Landsai-4 TM A MSS Data 

K. Yu 

Virginia Tech. 

An Analysis of the High Frequency Vibrations in Early TM Scenes 

J. Kogut 

Research & Data Systems 

Assessment of TM Band-to-Band Registration by the Block Correlation Method: 
Preliminary Results 

R.C. Wrigley 

NASA, Ames Research Center 




Landsat-4 Scientific Characterization 
Early Results Symposium 


TM Radiometry and Geometry Session, Cont. 
Thursday, 24 February 

Bldg. 8, Auditorium 
John Barker, Moderator 


3:05 p.m. Investigation of TM Band-to-Band Registration Using the JSC Registration Processor 

S.S. Yac 

Lockheed-EMSCO 

3:15 p.m. Geodetic Accuracy of Landsat-4 TM Data 

L. Oleson 

U.S. Geological Survey, EROS Data Center 
3:30 p.m. Adjourn 




Scientific Characterization 
Landsat-4 Early Results Symposium 

TM Applications Session 
Thursday, 24 February 

Bldg. J. Auditorium 
Darrel Williams, Moderator 


8:30 a.m. TM Quick-Look Analysis of the Washington, D.C. Area 

D. Williams 

Goddard Space Flight Center 

8:45 a.m. Preliminary Comparisons of Information Content Utility of TM vs MSS Data 

B.L. Markham 

Goddard Space Flight Center 

9:00 a.m. Information Content of Data from the Landsal-4 TM and MSS 

J. Price 

U.S. Department of Agriculture, ARS Hydrology Laboratoty 

9:15 a.m. Landsal-4 Investigation of TM and MSS Appliations 

D. Lauer 

EROS Data Center 
9:45 a.m. Break 

10:00 a.m. Preliminary Evaluation of TM Data for Renewable Resources 

D. Pitts 

NASA/Johnson Space Center 

10:20 a.m. TM Data Quality and Performance Assessment In Renewable 
Resources/ Agriculture /Remote Sensing 

R. Biuell 

Johnson Space Center 

10:50 a.m. Remote Sensing Investigations of Wetland Biomass and Productivity for Global 
Biosystems Research 

M. Hardisky 

College of Marine Studies, University of Delaware 

11:10 a.m. Assessment and Comparison of Landsat TM and MSS Data for Detecting Submerged 
Vegetation 

S. Ackleson 

College of Marine Studies, University of Delaware 


Landsat-4 Scientific Characterization 
Early Results Symposium 


11:30 a.m. 

11:90 a.m. 

12:10 p.m. 
1:30 p.m. 

1:90 p.m. 

2:10 p.m. 

2:30 p.m. 

*50 p.m. 

3:10 p.m. 


TM Applications Session, Continued 
Thursday, 24 February 

Bldg. 3, Auditorium 
Darrel William*, Moderate*' 

A Preliminary Evaluation of Thematic Mapper Data to Measure Suspended Sediments 
la Lake Chicot, Arkansas 

J. Ritchie 

U.S. Department of Agriculture 

Snow Reflectance from the TM 

J. Dozier 

University of California, Santa Barbara 

Lunch 

Preliminary Evaluation of T M :.oits 
D. Thompson 

NASA, Johnson Space Center 

Spacial Analysis of Thematic Mapper Products 

M. Jackson 

NERC, Experimental Cartography Unit, U.K. 

(Geologic Utility of Landsat-4 TM Data • Death Valley, California 

M. Ab»ams 

NASA, Jet Propulsion Laboratory 

Preliminary Study of Information Extraction of TM Data for a Suburban Regional 
Test Site 

D. Toll 

Goddard Space Flight Center 

Comparative Techniques Used to Evaluate TM Data for Land Cover Classification in 
Logan County, West Virginia 

H. Blodget 

Goddard Space Flight Center 

A Relative Accuracy Assessment of Landsat-4 MSS and TM Data for Land Cover 
Inventory 
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Y. C. Lu, 

Computer Science Corp. 



3:30 p.m. 


Adjourn 



Landsat-4 Scientific Characterization 
Early Results Symposium 

Fourth General Session 

Research Plans and Panel Discussion 
Thursday 24 February 

Bldg. S, Auditorium 
Mark Settle, Moderator 

This session makes available to the public a panel of Landsat-4 personnel to respond 
to questions on scientific, technical and policy topics. 


4:00 p.m. Thematic Mapper Research Plans 

Mark Settle 

Landsat-4 Program Scientist, NASA Headquarters 

4:20 p.m. Panel Discussion 

Panel members are: 

John Barker 

NASA, Goddard Space Flight Center 
Russ Koffler 

NOAA, National Environmental Satellite Information and Data Service 
Vince Salomonson 

NASA, Goddard Space Flight Center 

Mark Settle 
NASA Headquarters 


5:00 p.m. Adjourn 








LANDSAT-4 PROGkAM 
William Webb 

Landsat-4 Mission Operations Manager 
Goddard Space Flight Center 


The Landsat-4 program Is described In terms of its history, 
project objectives, mission requirements and system design. Accomplish 
ments, as well as future plans are outlined, Including the transition 
of the program to NOAA. 


NOAA's LANDSAT PROGRAM 


Russell Koffler 

Office of Satellite Data Processing 
and Distribution 

National Oceanic and Atmospheric Administration 


An Introduction to the National Oceanic and Atmospheric Administration's 
activities as the organization changed with the responsibility of 
managing the nation's operational civil remote sensing satellites. 
Emphasized will be NOAA's role In the Landsat Program. 

Discussion of related questions will take p^ce at the Symposium Wrap-Up 
Pai.el Session, 4:00 p.m., Thursday, 24 February. 


LANDSAT-4 SYSTEM DESCRIPTION 
T. C. Aepll 

Manager, Landsat Systems Engineering 
General Electric Co. 


The operation of the next generation Landsat system began following the 
launch of Landsat-4 on July 16, 1982. Involved are a new satellite de- 
sign, a new Thematic Mapper Instrument, new and varied data return paths 
and a new ground segment. The system has been operating for over six 
months durinq which time Initial assessments have been made of the perform- 
ance of these new elements. This paper summarizes the Landsat-4 system, 
highlights the new and different elements of Interest to data users and 
summarizes early results. 


INSTRUMENT DESCRIPTIONS 


Jack Engel 

Manager, Thematic Mapper Program 
Santa Barbara Research Center 
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LANOSAT-4 HIGHLIGHTS 

Darrel Williams 
Goddard Space Flight Center 

The first Landsat satellite was launched In July 1972. It carried an 
Earth-viewing sensor known as the Multlspectral Scanner (MSS) that obtained 
Imagery of the Earth's surface In four disc-els spectral bands. Landsat 1 
was followed by Landsats 2 and 3 In 1974 and 1978, respectively. Research 
conducted during the 1970's demonstrated that unique types of Information 
could be derived from MSS Imagery. For example, multlspectral Lar.dsat data 
has been used for geological mapping, crop forecasting, urban studies, land 
use planning, water management, map making, and a wide variety of other 
applications. More than one million Images of the Earth's surface have 
been collected by the first three Landsat satellites. This Image 
collection represents the first historical record of global surface 
conditions on our planet. Th< value of this Image collection has grown 
significantly during the first decade of Landsat oreratlons, and It will 
continue to grow In the future. 

Our ability to conduct routine observations of the Earth's surface on a 
repetitive, global basis has been considerably expanded with the successful 
launch of Landsat 4 In July 1982. Landsat 4 carries a new sensor know as 
the Thematic Mapper (TM) that Is far more sophisticated than the previous 
MSS Instruments. The TM measures the Intensity of surface radiation In 
seven discrete spectral channels (six of these channels are sensitive to 
variations In reflecced sunlight and the other channel, commonly referred 
to as the "thermal band". Is sensitive to surface temperature). Individual 
picture elements In a TM Image cover an area of 30x30 meters, except for 
the thermal channel which has a spatial resolution of 120x120 meters. 
Incoming radiation Is assigned to one of 256 gray levels In each spectral 
channel by the TM. In comparison, the MSS Is less sensitive with an 
overall range of 64 gray levels in each band. In summary, the TM possesses 
approximately twice the spectral resolution, three times the spatial 
resolution, and four times the sensitivity of the earlier MSS sensors. 
Considered together, the TM represents an order of magnitude Increase In 
our Earth observation capabilities from space. 


More than 100 Thematic Mapper Images have been fully processed and 
distributed to scientific Investigators. The preliminary results of 
current data analysis activities are described below. 

Agrlcul ture 

The TM spectral bands were selected In part for classifying crops and 
monitoring conditions. One of the earliest TM scenes was acquired over 
the Mississippi River In the central U.S. on August 22, 1982. A wide 
variety of crops are grown In this region, which encompasses portions of 
Arkansas, Missouri, Kentucky, and Tennessee. Analysis of this scene has 
demonstrated that several of the new spectral bands on the TM (particularly 
those at wavelengths of 1.6 and 2.2 micrometers) are particularly useful for 
distinguishing crops such as rice and soybeans, which were difficult to 
differentiate In earner MSS Imagery. Overall crop classification accuracy 
for this scene was Improved by roughly 25* using TM data. Early results 
suggest that a TM Image acquired at a single point in time may achieve the 
same accuracy In crop classification that has been achieved with multiple 
MSS images obtained throughout an entire growing season. The Improved 
spatial resolution of the TM Is particularly useful for crop assessment 
because a greater number of the TM picture elements fall within Individual 
fields. Picture elements (pixels) In single fields produce a more accurate 
measurement of the spectral characteristics of individual crops than do 
"mixed pixels" which straddle field boundaries and contain mixed spectral 
signatures of crops, roads, trees, water, etc. One of the early TM scenes 
acquired over the central plains of Canada has graphically demonstrated the 
utility of TM data for surveying areas that are cultivated In a "strip 
crop" fashion with alternating rectangular fields of cultivated and fallow 
land. Many portions of Europe, China, India, and South America contain 
agricultural fields that range In size from 10-40 acres. TM data should 
be particularly useful for crop surveys In these areas. 

Geology 

One of the other major factors In selecting the spectral channels on the 
Thematic Mapper was the desire of geologists to map variations in the 
abundance of clay minerals In semi-arid regions. Clay minerals generally 
possess diagnostic spectral absorption features at wavelengths of 2.0-2. 5 
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micrometers, and the TM band at 2.2 micrometers was specifically Included 
for this purpose. Geological evaluation of TM data has Initially 
concentrated upon a scene of Death Valley, California that was acquired 
November 17, 1982. Mlneraloglcal variations between major classes of rocks 
such as granites, quartzites, and volcanic lava:, can be clearly detected In 
the TM Imagery. Furthermore, subtle variations In clay mineralogy and 
abundance have been detected In the gravel fans that lie along the edges of 
the valley. These clay variations cannot be detected In aerial photography 
of the Death Valley area, demonstrating the utility of the new spectral 
channels at 2.2 and 1.6 micrometers for rock classification. A second 
major study site situated near Tucson, Arizona has also been examined An 
area of copper mineralization known as Silver Bell, Arizona has been 
studied In the past using airborne sensors that simulate the measurement 
capabilities of the TM. These TM simulator studies showed that certain 
types of clay minerals commonly associated with metal deposits could be 
spectrally detected. A TM scene acquired over the Tucson area by Landsat 4 
has confirmed the results of this earlier study, and an alteration "halo" 
of clay minerals associated with the Silver Bell deposit has been 
Identified In the orbital TM Imagery. Airborne TM simulator studies 
performed In the past were conducted over localized test sites roughly 
15x15 kilometers In size. The confirmation of these experimental results 
using early Landsat 4 data Implies that similar mineral evaluation studies 
can now be conducted routinely In semi-arid portions of the world, over 
areas of 185x185 kilometers or greater. 

Water Resources 

The Improved spatial resolution of the TM will lead to significant 
Improvements In our ability to map natural watersheds. It will Improve 
estimates of urban runoff and lead to better management of storm water In 
metropolitan areas as well. Snow and cloud-covered areas can be more 
accurately differentiated with TM data due to the unique spectral 
properties o f snow and clouds at wavelengths of 1.6 and 2.2 micrometers. 
These spectral channels will also provide a new capability for 
differentiating nutrients and sediments that are commonly foui'd In coastal 
waters. Preliminary analysis of TM data acquired over the southeastern 
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coastline of the U.S. near Charleston, South Carolina has demonstrated that 
suspended sediments can be distinguished from plankton In near-shore 
waters. Furthermore, the Improved spatial resolution of the TM should 
prove useful for longer term studies of topsoil erosion and sediment 
transport In coastal environments. The new thermal Infrared band on the TM 
will provide an Important new capability for monitoring soil moisture 
conditions and evaluating the effectiveness of crop Irrigation practices. 
Preliminary analysis has shown that surface temperature differences as 
small as 0.5° Centigrade can be detected In TM thermal Imagery. Large 
scale surveys of thermal conditions In vegetated regions will be useful for 
detecting water stress conditions In both natural and cultivated areas. 

Urban Studies/Land Use Planning 

The Thematic Mapper can be used to Identity surface features 30 meters on a 
side, which rougnly corresponds In size to a standard city lot. Highway 
construction, land excavation, urban growth, the health and extent of 
vegetation, etc. can be delineated In Thematic Mapper Imagery. Repetitive 
observations of these features Is an Important source of data for land use 
planning and resource management. The TM Is Ideally suited for collecting 
this type of data In a quantitative fashion at uniform scales and 
accuracies. Comparative analysis of TM Imagery acquired over large 
metropolitan areas such as Washington, DC and small rural towns such as 
Union City, Tennessee have demonstrated that a wide range of man-made 
features have unique spectral expressions In TM imagery. 

In summary. Initial analysis of TM imagery has produced highly encouraging 
results In a wide variety of disciplines. The TM possesses a unique set of 
measurement capabilities that will not be duplicated or surpassed by any 
orbital sensor system planned for the 1 980 ' s . Earth imagery acquired by 
the TM will extend and enlarge the value of the existing Landsat data 
collection. At the same time, the TM opens new avenues of Earth-related 
research that could not be addressed by earlier sensor systems. Evaluation 
and analysis of Landsat 4 TM data for NASA R&D purposes will continue 
during the next two years. The TM sensor will be transferred to NOAA in 
1985 for ooerational data collection. 


Dr. Mark Settle 

Landsat 4 Program Scientist 

NASA Headquarters 
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Dr. Vincent Salomonson 
Landsat 4 Project Scientist 
Goddard Space Flight Center 
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EARTH IMAGING AND DATA PROCESSING FOR MAPPING AND ANALYSIS, 
A TECHNOLOGICAL ASSESSMENT 


Ralph Bernstein 

IBM Palo A l4- o Scientific Center 
1530 Page Mill Road 
Palo Alto, CA 94303 0821 USA 


ABSTRACT 

Early observations of the earth used balloons and cameras. These platforms 
provided limited geographical coverage and small area images. The platforms 
were replaced by rockets and aircraft platforms, but did not provide repeat 
coverage or cost effective data. These Imaging systems used film, the 
procvdlng was chemical, and the Interpretation involved primarily human 
analysis of the photographic images. The current Landsat NASA/NOAA 
' *ndsat satellite program has the potential of providing users with data 
of exceptional spectral and dynamic range, and data of high resolution 
at nominal costs. In addition, the repeat coverage under controlled 
conditions provide data of significant scientific utility. Current sensor 
technology utilizes solid state detectors, and the data can be provided 
to users In a digital form. Thus, Information extraction and Interpretation 
can be aided by digital processing of the data. 

This paper describes early earth observation Imaging systems and compares 
the quality and utility of the previous data with the current Landsat-4 
Thematic Mapper data. In particular, the ability of the TM sensor to 
discern small natural and cultural features will be discussed and described 
with emphasis on land use analysis and mapping. The TM has a thermal 
band, capable of both day and night imaging. Both B&W and color products 
have been produced and will be presented and discussed. In addition, the 
wide spectral range of the TM data allow for new presentations of data 
to the users. Various band combinations In red/green/blue and Intensity/ 
hue/saturation will be presented, as well as experiments to attempt to 
display six bands in one presentation. 

Digital data processing was at one time the domain of data processing centers. 
This technology has not progressed to the point that Individuals now have 
their own computers. The paper will also describe an experiment that 
involves the use of an IBM Personal Computer to process the new Landsat-4 
Thematic Mapper data. This computer has been programmed to do statistical 
analysis of the data, produce graphical results, and perform limited image 
processing and display operations. This capability may be of importance 
to users who do not have access to larger computers. Future satellite 
communications technology, combined with the digital processing technology 
presents the possibility that future earth observation systems will provide 
data directly to the users, with a concurrent reduction in both time for 
data acquisition and data costs. 


INTERPRETARILITY OF CALIFORNIA'S AGRICULTURAL RESOURCES 
ON LANDS AT- 4 TM AND MSS IMAGERY--AN EARLY APPRAISAL" 


R. N. Colwell 
S. D. DeGloria 


Remote Sensing Research Program 
University of California 
Berkeley 94720 


ABSTRACT 


The purpose of this presentation is to provide an early appraisal of 
the interpretability of Landsat-4 TM and MSS image products for determining 
the extent to which major agricultural resources and conditions can be 
detected and identified. An evaluation of the quality of various spectral 
band combinations and transformations for a scene of the Central Valiev of 
California (Scene ID #84014518082X0, P: 43, R:34, 12/8/82) is accomplished 
through systematic image interpretation using photographic displays of the 
digital data, provided in part by IBM, in combination with low altitude 
oblique and ground photography and ground data provided in part by ongoina 
agricultural research projects in California under L'SDA sponsorship. Eva’- 
uation of the quality of the thermal band (TM6) for assessinq agricultural 
resources is also included in the interpretation process. The evaluation 
includes (1) a quantitative and qualitative interpretation of crop types, 
primarily small grains and fields und*»r double cropping regimes; (2) an 
assessment of the within-field, or feature, variability for assisting the 
interpretation process, and (3) an assessment of interpretation errors as 
influenced by image quality. 
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EVALUATION Of THEMAT 1C MAPPER PERFORMANCE 
AS APPLIEO TO HYDROCARBON EXPLORATION 
By Dr. John R. Everett, Dr. Charles Sheffield, Dr. Jon Dykstra 
Earth Satellite Corp., 1222 47th Street, Chevy Chase, MO 2061 5 


Results from a limited sample of Landsat-4 TM scenes demonstrate that 
the greater n-jmber of spectral bands and the Increased spectral and 
spatial resolution over previous satellite data greatly expand, Improve 
and refine the geological Inferences possible from space. The number and 
narrower width of available spectral bands, coupled with digital 
processing techniques such as band ratiolng, principal component analysis 
and hue-saturatlon-lntenslty presentation, allow the differentiation and 
mapping of a large number of geologic units and the Identification of 
altered rock associated with Intrusions In the Death Valley, California, 
area. These capabilities will contribute greatly to geologic exploration 
In arid areas. 

The Increased spatial resolution also assists In distinguishing 
glacial from tectonic features (difficult on previous Imagery) In the 
Detroit area. Several of the Ordovician oil and gas fields of Ontario lie 
on fracture zones clearly visible on Landsat-4 TM Imagery. All of these 
fields rely on fracture-related porosity and permeability tor production. 
The current-borne sediment plumes In Lake St. Clair and Erie demonstrate 
that Landsat-4 will contribute to our understanding of lacustrine, 
riverine and marine deposltlonal and sroslonal patterns and processes. 

Digitally enhanced Landsat-4 TM Imagery reveals geological features 
of the Anadarko Basin, Oklahoma, In greater numbers and detail than 
previously possible. Comparison with existing work suggests that some of 
the subtle spectral differences detected are related either to changes In 
chemical composition or resultant changes In vegetation related to 
micro-seepage of hydrocarbons adjacent to oil fields In the area. 

There seem to be several categories of features (throughgol ng 
structures, large area subtle tonal contrasts, etc.' that are still best 
seen on older Landsat Imagery. Thus, Landsat-4 does not make earlier 
systems obsolete; rather It provides a powerful new tool to exploration 
geo log I sts. 


GEOLOGIC UTILITY OF LANDSAT-4 TM DATA: 
DEATH VALLEY, CALIFORNIA 

Michael Abrams, Anne Kahle, Alan Gillespie, 
James Conel , Harold Lang 
Jet Propulsion Laboratory 
Pasadena, CA 


The present LIDQA study Is being conducted In several phases to quantify 
the performance of the TM *»i s -a-v1s various geological applications. These 
phases Include: (1) analyses of the geological utility of the data with 

respect to the Increased spatial resolution and number of bands (compared 
to the MSS); (2) analysis of geometric accuracy; (3) analysis of radiometric 
performance of the TM scanner. 

This report presents some preliminary results of geological analysis of 
TM data for Death Valley, California, and observations on TM data quality. 
Geometrically and radlometrlcal ly corrected data (CCT-PT) have been 
examined for scene E-40124- 1 7495 , acquired November 17, 1982. Overall, the 
TM data appeared to contain a marked Increase In geologically useful 
Information; however, a number of Instrumental or processing artifacts 
may well limit the ability of the geologist to fully extract this Information. 

A 27 by 18 km area was extracted from the scene for detailed analysis. 
Color-enhanced band composites, principal components composites, and color 
radio composites were produced from the TM data. Lithologic Interpretation 
maps were prepared and compared to detailed published geologic maps. More 
detail and delineation of alluvial fan units was evident on the TM data; 
most of the sal t/sl 1 t/evaporl te units In the Death Valley playa were 
apparent on the Images; many of the outcropping volcanic and sedimentary 
rock types were separable based on differences in their spectral response. 
However, substantial shadowing due to the combination of low sun angle and 
severe topographic relief limited the ability to Intepret information In 
the mountainous areas. High sun angle, summer scenes would greatly reduce 
this problem. 

The same region was extracted from Landsat-2 MSS data and also from aircraft 
mul tlspectral thermol data. Displayed as color composites and principal 
component Images, the TM data were superior to the MSS data In all respects. 
The multlspectral thermal data, however, allowed separation of alluvial 
units and lithologic units not seen on the TM data. This indicates the 
complementary nature of the two data sets for providing lithologic 
separations based on sensing different physical properties. 

Two significant radiometric problems were observed In the TM data. The 
mean signal levels differ appreciably between the forward and reverse 
scan directions, possibly due to inaccuracies In the radiometric calibrations. 
This effect manifests itself as banding in harshly contrast-enhanced images. 
The second problem is most pronounced in ratio images and appears as 
striping, most severe at the edges of the scenes. The striping consists 
of paired lines of high-variance noise, with a 17-line periodicity. They 


occur at the seams between the forward and reverse scan swaths, where 
extension scan lines have been Interpolated to fill data gaps. It appears 
that the Interpolation procedure used to process this scene was faulty. 

The TM data were registered to a topographic base map. Results of the 
registration procedure Indicate that the error In geometric accuracy of 
the CCT-PT data Is less than 2 pixels. A more rigorous test of accuracy 
Is planned. 
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AN INITIAL ANALYSIS OF LANDSAT-4 THEMATIC MAPPER DATA 
FOR THE DISCRIMINATION OF 

AGRICULTURAL, FORESTED WETLAND, AND URBAN LAND COVERS 

By: Dale A. Quattrochi 

NASA/Earth Resources Laboratory 

Code HA30 

National Space Technology Laboratories 

NSTL Station, MS 39529 

ABSTRACT 

An Initial analysis of Landsat-4 Thematic Mapper (TM) data for the 
discrimination of agricultural, forested wetland, and urban land covers 
was conducted using a scene of data collected over Arkansas and Tennessee 
on August 22, 1982. A study area in Poinsett County, Arkansas was used 
to evaluate a classification of agricultural lands derived from multi- 
temporal Landsat Multispectral Scanner (MSS) data in comparison with a 
classification of TM data for the same area. Results indicated that th. 

TM performed significantly better than the MSS in cverall classification 
accuracy. Data over Reel foot Lake in northwestern Tennessee were used to 
evaluate the TM's capabilities for delineating forested wetland species. 

A classification derived for the study area yielded an overall accuracy of 
95% correct in discriminating five categories of forested wetlands. Addition 
ally, the TM data were used to identify urban features within a small city 
in northwest Tennessee. A classification of Union City, Tennessee was evalu- 
ated for accuracy in discriminating urban land covers. This classification 
produced an overall accuracy of approximately 90% correct when compared with 
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ground truth Information for the study area. An assessment of dlqltall.v 
enhanced TM data was also performed using Drlnclpal components analysis 
to facilitate photointerpretation of urban features of the Union City area. 
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Canadian Plans for Thematic Mapper Data 


Authors W.M. Strome, F.E. Guertln, A.B. Collins and D.G. Goodenough 

Presented by W.M. Strome 
Canada Centre for Remote Sensing 
2464 Sheffield Road 
Ottawa, Ontario, Canada 
K1A 0Y7 


ABSTRACT 


Canada was the first country outside the U.S.A. to receive, 
process and distribute data from the original LANOSAT series. The data 
from these satellites have been used to Improve the state-of-the-art of 
resource management and environmental monitoring In Canada over the 
past decade. The utility of MSS data has been proven In a wide range 
of applications, from Ice reconnaissance, through agriculture and 
forestry to geological exploration, to name a few. Despite the 
usefulness of MSS, Its limitations prevented full use of satellite 
remote sensing technology for some purposes. Thus, It was decided to 
upgrade the Canadian facilities to receive and process the data from 
LANDSAT-4 . 

The modifications to the existing equipment for processing MSS 
data from the new satellite were relatively straightforward and are now 
complete. The Prince Albert Satellite Station has been modified to 
permit reception and recording of TM data for weitern North America. A 
small scale system has been built to transcribe TM data from the high 
density station tapes to computer compatible tape. The throughput of 
this system Is very low - two to five scenes per week. However, It 
does enable us to develop and test the algorithms to be used In our 
Thematic Mapper Bulk Processing system and our Multi Observation 
Satellite Image Correction System (MOSAICS). These will be operational 
in 1984 and 1986 respectively. The latter system will provide geocoded 
Image data, keyed to maps and Independent of satellite orbit and sensor 
resolution. To analyze the TM imagery a new facility called the 
LAHDSAT-D Image Analysis System Is being developed and will be phased 
gradually Into operation starting In 1983. 


LANDSAT-4 IMAGE DATA ANALYSIS 


Robert B. MacDonald 
NASA/Johnson Space Center 
Houston, Texas 


Researchers at the Johnson Space Center have completed a preliminary 
analysis of data acquired by the Thematic Mapper over several test sites. 
These data were analyzed, together with "site ground truth measurements" 
to evaluate various performance characteristics of the Thematic Mapper. 
The results of this analysis are presented. 


Third General Session 

Scientific Characterization 
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LANDSAT-4 EVALUATION PROGRAM AND SCIENTIFIC 
CHARACTERIZATION ACTIVITIES 


John L. Barker 

Landsat-4 s* sod ate Project Scientist 
Goddar*' Space Flight Center 

As part of the Landsat-D program, NASA committed to the following 
objectives: 

1) Thematic Mapper (TM) Capability Assessment 

2) Multi spectral Scanner (MSS) to TM Transition 

3) Operational System Feasibility Demonstration 

4) Continuity of MSS Imagery 

5) Continued Foreign Access 

With the turnover of Landsat-4 MSS operations to NOAA on January 31, 
1983, NASA has effectively accomplished 3 of the 4 engineering objectives 
namely 3, 4 and 5. Objective 2 Is being addressed with the development 
of the TM Image Processing System (TIPS). Objective 1 falls primarily 
Into the pervlew of scientific characterization rather than engineering 
verification of reliability, specifications and quality assurance. The 
characterization objectives of the Landsat-4 Science Office at Goddard 
Space Flight Center are: 

1) Characterize Accuracy and Precision of Sensor and Spacecraft 
Performance 

2) Characterize Accuracy and Precision of Image Data Quality 

3) Characterize Accuracy and Precision of Derived Information 

4) Recommend Landsat-4 System Improvements 

5) Communicate Results to Research Community 

In order to accomplish these objectives, NASA has embarked on a three 
pronged effort. Firstly, In-house activities at Goddard have been 
directed toward full access and utilization of the pre-launch and 
In-orbit engineering test data on the sensor and spacecraft. Secondly, 


the principle scientists In land remote sensing have been gathered together 
as part of a major scientific characterization effort following the LIDQA 
Application Notice (AN) In 1981. Two Landsat-4 Workshops have already 
been held for these Investigation teams, namely In May, 1981 and 
January, 1982. Most of the Investigators are now under contract. Investi- 
gators who have already received Imagery are presenting at this Landsat-4 
Early Results Symposium. A final symposium Is scheduled for the Fall of 
1984, prior to the turnover of TM to NOAA. In a third phase, NASA Intends 
to move from an emphasis on sensor and spacecraft performance and image 
data quality Into a more extensive scientific Investigations program. 
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MSS EUiUAT.ON REPORT 


William Alford 

landsat.-4 Senior Systems Programmer 
Goddard Space Flight Center 
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TM REPORT 


John L. Barker 

Landsat-4 Associate Project Scientist 
Goddard Space Flight Center 

In this TM summary report, emphasis Is placed on the key spectral, 
radiometric and geometric characteristics of the sensor. Future editions 
will emerge from the studies of the Landsat-4 LIDQA Investigations team 
and will result In an expanded emphasis on image data quality and appli- 
cations Information. 

Spectrally the key findings of Importance to scientific use of TM 
In the visible (VIS) and near Infrared (NIR) Included narrower bands, 
and new bands in the shortwave and thermal Infrared (SI*IR and TIR). 

Band 5 Is slightly wider than specified and Its upper edge at 50% of 
maximum of 1784 nm (nanometers) includes some unintended overlap with 
water absorption bands. The thermal Band 6 is half the expected width, 
which causes no foreseeable scientific problems, especially since it 
has nearly twice the expected signal-to-nolse ratio. 

Radiometrical ly, the TM sensor Is characterized by extreme linear- 
ity and stability in comparison to previous MSS photomultiplier systems. 

The forward and reverse scans on TM have shown some unexplained radio- 
metric differences of up to four digital levels, especially in the 
four bands on the primary focal plane ( PFP ) . Intraband variability 
between detectors has been less than about ± 0.3% for the four PFP 
bands and about ± 1% for Bands 5 and 7. This means that radiometric 
calibration of Bands 5 and 7 is more important for within image de- 
striping. PFP bands, as monitored by the internal calibration (IC) have 
shown a monotonic decrease in gain of about 5% in the five months since 
launch. On the other hand, the relative internal gains of Bands 5 and 7 
have varied with a period of about 60 days and an amplitude of about 6%. 
This difference in pattern between PFP and the cold focal plane (CFP) bands 
suggests a band-dependent gain change that is independent of possible IC 
changes with time. This indicates that it is essential to apply IC radio- 
metric calibration before any histogram equalization algorithms in order to 


preserve both relative and absolute Interband radiometric accuracy. 
Slgnal-to-nolse ratios exceed specifications by nearly a factor of two 
In all bands. Observed background noise Is ± 1 count In Bands 1, 5 and 
7 and about t .5 counts In the other three bands. Detector 3 In Band 5 
Is dead. In general, the first detectors are the noisiest. Channel 7 
In Band 7 Is more than a factor of two noisier than other detectors, 
and hns shown a 20% gain Increase since pre-launch calibration In 
system-level thermal vacuum testing. Thermal Band 6 showed a loss In 
gain of approximately 40% from Augusc to January, when the CFP cooler was 
outgased and Band 6 regained Its Initial In-orbit gain sensitivity. 

This gain loss corresponds approximately to a change of from 0.3 to 
0.45°C per digital count. There appears to be a mislabeling of the Band 
6 channels for the calibration data. Assuming that there Is no signifi- 
cant limitation from quantization, dynamic ranges may be too large for 
typical scenes, especially In Bands 6 and 7. 

Geometrically, pre-launch tests indicated compliance with most 
specifications. Channel 4 In Band 2 has a smaller modulation transfer 
function (MTF ) than desired and Is being replaced by Its neighbor prior 
to geometric resampling. There Is an apparent PFP-CFP misregistration 
of .75 pixels along scan and 0.2 pixels across-scan. An Intermittant 
geometric discontinuity of up to 2 pixels between forward and reverse 
sweeps has been observed and Is under investigation. Peak-to-peak 
attitude deviations are being measured by the angular displacement sensor 
(ADS) at ;s than 0.5 pixels and are removed during systematic geometric 
preprocessing. 



ANALYSIS OF THEMATIC MAPPER DATA FROM CEOLOCIC 
TERRAIN ’HITS IN A CALIFORNIA-NEVADA SCENE 


Rupert Haydn 

Zentralstelle fur Geo-Photogrammetrie w. Fernerkundunq 
Luisenstrasse 37, 8000 Munchen 2 
Federal Republic of Germany 

and 


Nicholas M. Short 
Geophysics Branch 
NASA/Goddard Space Flight Center 
Greenbelt, MD 20771 


Abstract 


Thematic Mapper data from the Death Valley scene (Landsat 40124- 
17495; November 17, 1982) have been processed for a variety of output 
products. Edge-enhanced natural and false color composites of the full 
and partial scene(s), using the intensity, hue, saturation (IHS) color 
system, were produced on a DEC PDP-11. 

Ratios of selected reflectance bands including the Short Wave 
Infrared bands 5 and 7, can be displayed in combinations in which two 
or three ratios are expressed in the IHS color scheme and another ratio 
is used to calculate pixel shifts by which an image pair is created to 
generate a 3-D stereo effect. This approach can also depict variations 
in temperature, from band 6 data, as the stereo component. Examples of 
these experimental products will be displayed in 3-D to the audience. 

A data subset from the northern Death Valley area was analyzed 
on Goddard's IDIMS computer. Unsupervised classifications, based on 
a maximum likelihood classifier, yielded nine separable classes from TM 
bands 2, 3 and 4 (roughly equivalent to the MSS interval) and 11 classes 
for bands 1, 2, 3, 4, 5 and 7 combined. Discrimination of rock types 
in the mountain terrain was poor but recognition of discrete units 
in the valley was notably improved. A supervised classification 
based on twelve defined classes including five rock units present in the 
mountains was superior to the unsupervised ones. Better definition 
of valley units was achieved with color composites made from several 
image combinations of principal components derived from the Karhuuen- 
Loeve transformation function developed by ERIM to specify brightness 
and greenness. 
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TM SENSOR RADIOMETRY 
John L. Barker 

Landsat-4 Associate Project Scientist 
Goddard Space Flight Center 

An overview of the pre-launch and In-orbit radiometric 
calibration procedures Is provided. Including both absolute and 
relative calibration. 



SENSOR GEOMETRY 


J. Engel 

Santa Barbara Research Center 



TM RADIOMETRIC PROCESSING 


John L. Barker 

Landsat-4 Associate Project Scientist 
Goddard Space Flight Center 


An overview of the SCROUNGE-Era TM radiometric digital processing 
Is provided from abstraction of calibration data In ADDS througn radio- 
metric calibration of the CCT-BT to CCT-AT tapes and subsequently 
through geometric resampling to CCT-PT tapes. 


TM INTRODUCTORY SESSION 


GEOMETRIC PROCESSING 
Eric Beyer 
General Electric Co. 


An overview of TM Geometric Processing Is provided. The description of 
the Flight Segment Includes the TM, attitude measurement sensors, attitude 
control system and payload correction data. The ground processing description 
includes payload correction processing (In which Flight Segment telemetry Is 
converted Into geometric distortion Information), ground control point pro- 
cessing and Image data reformatting and resampling. 



IMAGE DATA QUALITY 


J. Lyon 

Goddaro Space Flight Center 


MSS Session 





MSS RADIOMETRIC ACCURACY ASSESSMENT 


William Alford 

Landsat-4 Senior Systems Programmer 
Goddard Space Flight Center 




Spectral Characterization of the Landsal-4 MSS Sensors 


ORIGINAL RAU: IS 
OF POOR QUALITY 


Brian L. Markham and John L. Barker 
NASA/Goddard Space Flight Center 
Earth Resources Branch 
Greenbel t, MO 20771 


ABSTRACT 

Relative spectral response data for the Landsat-4 and landsat-4 backup multispectral 
scanner subsystems (MSS), the protoflight and flight models, respectively, 
are presented and compared to similar data for the Landsat 1, 2, 3 scanners. 
Channel-by-channel (six channels per band) outputs for soil and soybean targets 
were simulated and compared within each band and between scanners. The two 
Landsat-4 scanners proved to be nearly identical in mean spectraj response, 
but they exhibited some differences from the previous MSS's. Principal differences 
between the spectral responses of the Landsat-4 scanners and previous scanners 
were: (1) a mean upper-band edge in the green band of 606 nm compared to 
previous means of 593 to 598 nm, (2) an average upper-band edge of 697 nm in 
the red band compared to previous averages of 701 to 710 nm, and (3) an average 
bandpass for the first near-IR band of 702-814 nm compared to a range of 693- 
793 to 697-802 nm for previous scanners. These differences caused the simulated 
Landsat-4 scanner outputs to be 3 to 1C percent lower in the red band and 3 
to 11 percent higher in the first near-IR band than previous scanners for the 
soybeans target. Otherwise, outputs from soil and soybean targets were only 
slightly affected. The Landsat-4 scanners were generally more uniform from 
channel to channel within bands than previous scanners. One notable case of 
poor uniformity was the upper-band edge of the red band of the protoflight 
scanner, where one channel was markedly different (12 nm) from the rest. For 
a soybeans target, this nonuniformity resulted in a within-band difference of 
6.2 percent in simulated outputs between channels. 
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INVESTIGATION OF RADIOMETRIC PROPERTIES OF LANDSAT-4 MSS 

Daniel P. Rice 

Environmental Research Institute of Michigan 
Ann Arbor, Michigan 48107 


Summary 


O bjectives 

The objectives of our investigation a r e to characterize Landsat-4 
multispectral scanner (MSS) image data quality relative to: 

(a) Detector calibration: Study the calibration of the six detec- 

tors in each band in order to determine the magnitude of any calibration 
differences that remain after ground processing, and, If needed, provide 
Information that would support corrective techniques. 

and (b) Satel 1 ite-to-satel 1 Ite calibration: Study calibration dif- 

ferences between Landsat-4 and previous Landsat satellites and, as needed, 
develop a metnod to adjust Landsat-4 multispectral scanner (MSS) signals, 
in all four spectral bands, to match the calibration of previous MSS 
sensors. 


General Approach 

Two full frames of radiometrical ly corrected (Type A) Landsat-4 MSS 
data were analyzed digitally and visually for evidence of residual cali- 
bration differences between detectors, quantization effects, and other 
sensor-related artifacts. Both standard statistical and Fourier analysis 
techniques were employed. 

Opportunities for coincident coverage by Landsats 3 and 4 were 
identified in the contiguous 48 states. Paired acquisitions for two 
scenes have been obtained to date and analyzed to establish relation- 
ships between signal values from common areas imaged by the two MSS 
systems. 
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I nitial Findings 

General Radiometry . We found that the Landsat-4 MSS produces data 
of generally good quanty with dynamic ranqes and target responses simi- 
lar to those of previous Landsat systems (See Figure 1). Clouds often 
caused saturation in all bands. The extent of specific system artifacts 
are discussed individually below. 

Detector-to-Detector Differences . Striping caused by residual 
differences in detector cal i brat ion and quantization effects appears 
to be quite well corrected. The histogram equalization algorithm used 
to adjust radiometric look-up tables appears to be working as intended. 
Residual rms differences of 0.3 digital counts or less were measured 
between detector means from a variety of scene areas. Fourier analysis 
of a down-track profile obtained by averaging the pixels in each scan 
line showed amplitudes attributable to the residual detector differ- 
ences that were also approximately 0.3 digital counts or less. These 
differences were indicated by disturbances in the transforms at wave- 
lengths of two, three, and six scan lines. 

Coherent Noise . A low-level coherent noise effect was found in all 
bands, appearing in uniform areas as a diagonal striping pattern. The 
principal component of this noise was found by Fourier analysis to be at 
a consistent wavelength of 3.6 pixels along the scan line (See Figure 2), 
corresponding to a noise frequency of about 28 KHz. The magnitude varied 
from about 0.75 of one count in the worst case (Band 1) to 0.25 counts 
in the best case (Band 4). This noise pattern was visually discernable 
in water areas of a Landsat-4 image but not discernable in a simultaneous 
Landsat-3 image of the same area. 

Landsat-4 to Landsat-3 Calibration . Initial regression relation- 
ships were computed between' average signal amplitudes from a number of 
areas in common between Landsat-3 and -4 images of the two scenes with 
paired acquisitions. Linear fits with R? values exceeding 0.99 were 
obtained (See example in Figure 3). Letting 

Landsat 3 = Gain x (Landsat 4) + Offset 

we computed gains that were higher (by approximately 10% in Bands 1, 2, 
and 4) than ones based on prelaunch calibration values for the two 
sensors. While 23 points were used to determine this initial result, 
updated relationships based on more data points will be presented in 
the complete paper. 
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Evaluation of Landsat-4 Multi spectral 
Scanner Ground Segment Performance* 

N. A. Bryant, A. L. Zobrist, F. C. Billingsley, 

S. Z. Friedman, B. Gokhman, T. L. Logan 

An analysis was performed on the P-Format computer compatible tapes for the 
Landsat-4 Multlspectral Scanner (MSS) scene 40109-15140 acquired November 2, 
1982 over Washington, D. C. and environs. Three tests of sensor geometry were 
undertaken: a) band-to-band registration, b) llne-to-line registration 
between swaths, and c) geometric correction of sensor and spacecraf t/ephemerls 
characteristics. The band-to-band and llne-to-line registration was measured 
at one hundred pixel spaclngs along a line using the phase correlation Image 
alignment method developed by Kuglin and Hines and adapted to a one- 
dimensional FFT correlation technique. The sensor and spacecraft geometric 
calibration analysis was checked by registering ground control point chips 
developed under the Landsat-3 Master Data Processor program. Twenty four 
ground control points (gep) were located In the scene using the phase 
correlation Image alignment method. (Note: this contrasts with 2 geps being 
found by the Landsat-4 ground data processing segment using an autocorrelation 
technique.) A least squares fit of the geps was computed and the vector 
offsets of the residuals was tabulated and plotted. 

An analysis of the tabulated and plotted results revealed the following 
characteristics: a) There is no measurable swath-to-swath misregistration in 
the P-format Imagery, b) Band-to-band registration is 0.1 pixel or less on 
the average. c) Line-to-line jitter effects may be occurring, but could not 
be more than one half pixel offset In any instance, d) It is apparent from 
the line-to-line plots that approximately 1/8 pixel shift to the right 
predominates. This may be a scene dependent artifact associated with the 
eleven degree off-polar orbit, an artifact of our phase correlation algorithm, 
or indicate a need to tune the estimate of spacecraft velocity for the 
geometric calibration. e) There appears to be a brightness modulation with a 
pseudo coherent pattern that gives a false Impression of swath-to-swath 
misregistration. The brightness modulation is more severe on the right than 
left side of the frame, and may be associated with Jitter impacting dwell time 
rates over pixels. f) Ground control point offsets from the residual surface 
fit are generally good ( RMSE of 2.5 pixels or 142.5 meters), although some 
points are off to a significant degree. An analysis of the TM P-data revealed 
a slow spacecraft roll condition under one hertz in this scene. As only two 
geps were found by the Landsat-4 ground processing system, this systematic 
distortion was not adequately removed. 

•This abstract presents the results of one phase of research carried out at 
the Jet Propulsion Laboratory, California Institute of Technology, under 
contract No. NAS7-100, sponsored by the National Aeronautics and Space 
Administration. 
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RADIOMETRIC CALIBRATION AND GEOCODED PRECISION PROCESSING 
OF LANDSAT 4 MULT I S PECTR AL SCANNER PRODUCTS 
BY THE CANADA CENTRE FOR REMOTE SENSING 


AUTHORS: J. MURPHY, D. BENNETT, F. GUERTIN 


PRESENTED BY JENNIFER MURPHY 


CANADA CENTRE FOR REMOTE SENSING 
2464 SHEFFIELD ROAD 
OTTAWA, ONTARIO 
CANADA 


ABSTRACT 

The method used by the Canada Centre for Remote Sensing 
(CCRS) for the radiometric calibration of Landsat 4 
M u 1 1 i s pe c t r a 1 Scanner (MSS) data is reviewed, with reference 
to the methodology used by the Centre for Landsat 1 2 and 3 
MSS data Inherent in this technique is the possibility for 
the user to convert the corrected digital values to the 
absolute scene radiance of the target under observation. 
The generation of the constants needed for this final 
conversion requires both tne pre-launch and post-launch 
radiometric calibration constants as supplied by NASA. 
Results of some preliminary comparative studies of the 
radiometric properties of the Landsat 4 MSS versus earlier 
satellites in the Landsat series are presented. In 
addition early observations on the stability of the 
calibration data firstly within one scene, secondly, 
within one orbit, and thirdly, over a period of several 
months, are presented. Quantitative estimates of residual 
striping in the corrected products are also presented. 

The method used by CCRS, to perform precision processing of 
Landsat MSS data to generate geocoded or cap compatible 
Landsat MSS products in the Universal Transverse Mercator 
projection, is reviewed. Landsat 4 MSS precision processed 
products are evaluated for geodetic accuracy, and are 
compared to similar products from the previous Landsat 
satellites to assess the orbit independent registration 
accuracy . 
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Landsat-4 Mult Ispectral Scanner Radiometric Intraband Performance 
Eric Ellason, Hugh Kieffer & Kay Edwards 
U. S. Geologic Survey, Flagstaff, Arizona 86001 

A geometrically raw Image of Washington D.C. was aqulred from EROS Data 
Center. These data have been radlometrlcal ly corrected by GSFC; radlo- 
metrically raw data are not available. The data aqulred show little of the 
detector striping common In earlier Multlspectral Scanner Images. 

The radiometrlcal ly corrected data have uniform means and standard 
deviations for the detectors In each band; however, the data for different 
detectors utilize a different pattern of DN levels, (typically 1 DN level out 
of 3 Is not used), resulting In ubiquitous striping of 1 DN amplitude. 

All bands have pronounced coherent noise obvious In contrast enhanced 
Images, and evident In natural color composites. One-dimensional Fourier 
transforms of a flat field show a sharp peak with a period of 3.5 samples. 
There Is evidence of weaker period noise at other frequencies. We have 
attempted to remove the coherent noise by both transform techniques and by a 
specially designed sequence of "box" filters. 
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LANDSAT-4 SCENE-TO-SCENE 
REGISTRATION ACCURACY ASSESSMENT 

f f'AL PAGE IS 

TY By: James E. Anderson 

NASA/Earth Earth Resources Laboratory 
Code HA20, Bldg 1100 
NSTL Station, MS 39529 

ABSTRACT 

This report documents initial results obtained frcm the registration of 
Landsat-4 MSS data to Landsat-2 MSS data. A comparison is made with results 
obtained from a Landsat-2 MSS-to-Landsat-2 MSS scene-to-scene registration 
(using the same Landsat-2 MSS data as the "base" data set in both procedures). 
RMS errors calculated on the control points used in the establishment of 
scene-to-scene mapping equations are compared to errors computed from inde- 
pendently chosen verification points. Models developed to estimate actual 
scene-to-scene registration accuracy based on the use of electrostatic plots 
are also presented. This project will include analyses of TM data at a later 
date, and both SCROUNGE and TIPS era products will be evaluated. 
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GEOMETRIC ACCURACY OF LANDSAT-4 IMAGE DATA 

R. Welch and E. Lynn Usery 
Department of Geography 
University of Georgia 
Athens, GA 30602 


Landsat-4 MSS digital data received from the EROS Data Center (EDC) 

In CCT-P formats are being rectified with the aid of a modified version 
of the Digital Image Rectification System (DIRS) originally developed at 
the NASA Goddard Spaceflight Center. The rectification procedure involves 
the digitization of ground control points (GCP's) from U.S. Geological 
Survey 1:24,000 scale topographic maps and the determination of Image 
locations with the aid of an interactive image processing system. Once 
the coordinates of the GCP's are determined In both the map (UTM) and image 
(pixel and line) coordinate systems, these data ar*» transferred to an 
IBM 370/158 mainframe computer on which the DIRS software Is resident. 
Corrections based on a least-squares squares fit o f polynomials of the 2nd 
to the 5th degree are developed and applied to the Imago coordinates. 

Accuracy checks are performed by two methods: 1) point pair distance 

checks involving the comparison of distances computed from the map and 
image coordinates; and 2) vector plots of the errors of withheld GCP's to 
which the least squares polynomial equations have been applied. 

Results of analyses of an MSS scene of a rugged area in North Georgia 
have provided the following results. Distances between point pairs of 
GCP's computed from the (raw) CCT-P data provided by EDC produced errors 
of approximately 120-130 m. Thi* indicates a relative positional error of 
about 50-60 m at well-defined points. After rectifying the full scene 
using 42 GCP's and 5th degree polynomial s, a root-mean square error (RMSE) 
of +40 m in planimetric position was obtained. Subscene areas of single 
1:24,000 scale quadrangles which have been rectified produce errors of 
about 30 m between point pairs. 

The procedures for locating the GCP's with the interactive image 
processing system are bei lg revised to allow sub-pixel determination of 
image coordinates. This refinement should yield higher accuracies. Work 
to-date indicates that the largest error source is the difficulty in locating 
the GCP's in the MSS scenes, and not defects in image geometry. 
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GEODETIC ACCURACY OF LANDSAT-4 MULTI SPECTRAL SCANNER AND THEMATIC MAPPER DATA 
Principal Investigator: June M. Thoroodsgard 

EROS Data Center, U.S. Geological Survey 


The EROS Data Center, as the primary facility for generation and dissemination 
of Landsat data products, is evaluating the geodetic accuracy of data from both 
the Multlspectral Scanner (MSS) and Thematic Mapper (TM) processing systems. 

The standard Landsat-4 MSS and TM digital tapes contain information which 
relates the map projection coordinate system to the image coordinate system. 
Therefore, by converting the latitude and longitude of a given location to the 
map projection coordinates, the image coordinate location can be calculated. 
Verification of the Image's geodetic accuracy can then be accomplished by 
comparing the calculated image coordinates with the visually determined image 
location. 

Landsat-4 digital images produced by NASA are either system-corrected or 
precision-corrected (ground control points applied). This method of geodetic 
verification has been applied to images of both precision-corrected and 
system-corrected MSS and TM data. 


MSS GEOMETRIC ACCURACY ASSESSMENT 


M. Imhoff 

Goddard Space Flight Center 


MSS GEOMETRIC PERFORMANCE: 

L I NE-TO-L INE DISPLACEMENT ANALYSIS 


L. Fusco 

Earthnet Program Office 
Frascati, Italy 
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Impact of Landsat MSS Sensor Differences 
On Change Detection Analysis 


ft* 


William Likens and Robert Wrigley 
NASA Ames Research Center 
Moffett Field, California 94035 


INTRODUCTION 

This work examines MSS sensor radiometric and geometric 
differences for their impact upon change dfection. Comparative 
response of various MSS sensors becomes a consideration in change 
detection analysis. Differences in overall sensor response, due 
to both factors within the sensor system, and external factors 
such as atmospheric effects, must be dealt with before change 
detection techniques using image differencing can be successfully 
implemented . 

To date three MSS scenes have been acquired and compared for the 
San Francisco area: 1) October 12, 1982 Landsat 4, 2) October 8, 
’982 Landsat 3, and 3) October 4, 1981 Landsat 2 data. These 
data comprise a common set that has been intercompared . Three 
additional comparison sets have been ordered and will scon be 
examined: 1) Landsat 2 and 4 data of Artesia, New Mexico, both 
acquired on October 9, 1982, 2) Landsat 3 and 4 data of Rhode 
Island, both acquired on December 22, 1982, and 3) Landsat 3 and 
4 data of New Hampshire, both acquired on December 22, 1982. All 
the data are being acquired in the EDIPS A-Tape (radiometric 
corrections only) format. While verbal communication with 
Goddard and EROS staff had given the impression that data for the 
same area would not be acquired concurrently from Landsat 
satellites, a search of the EROS Landsat data base showed that of 
the approximately one dozen possible instances of concurrent 
coverage over the United States between September and December, 
two instances existed where concurrent coverage was aquired. The 
imagery ordered were those scenes having the least cloud cover 
for the two orbit paths for which concurrent coverage existed. 
The most useful results of this investigation are expected from 
examination of the concurrently acquired image pairs, rather that 
from examination of the San Francisco data. 

ANALYSIS TO DATE 

The method of MSS comparison being used here is to coregister the 
data to a base, followed by between image scattergram generation 
and evaluation, and image differencing. In this case the Landsat 
2 and 3 data have been registered to the Landsat 4 data. The 
area coregistered was limited to a small portion of the scene 



near the city of Santa Cruz because this was the on±y area cloud 
and haze free on all the images. Visual examination of each 

image, as well as of the between image scattergrams , are intended 
to provide information of the radiometric characteristics of each 
sensor. The image differencing, after normalizing for 
radiometric differences, provides a visual depiction of geometric 
distortions . 

Comparison of the San Francisco data sets show that Landsat 4 
appears to have no problems that would preclude comparison to 
other MSS data for change detection, though there are comparison 
considerations that are common to comparison of any MSS data that 
must be considered. The Landsat 4 data was noted t.o have 
diagonal banding in bands 1, 2, and 3. This type of banding, or 
"woodgrain" effect has been noted to occasionally occur with the 
other MSS sensors. A subjective conclusion was that in this 
instance the magnitude of the effect was not great enough to pose 
a significant problem, though this requires and will receive 
further evaluation when additional image pairs are examined. 
Empty bins, as a result of the signal decompress.* r. ; , process, were 
noted in Band 4, though NOAA now reports this has now been 
corrected. Malila and Rice have previously reported both the 
diagonal banding and the empty bin effects at the LIDQA meeting 
held in January 1983 at GSFC. Band 2 of the Landsat 3 data was 
misregistered with respect to the other bands in the bottom 
portion of the scene. We believe this to be a ground processing 
rather than a sens-or problem, and have noted this effect once 
before in Landsat 3 EDIPS processed data set (described then as a 
"dogleg"; both a stretch and a leftward shift increasing towards 
the bottom of the image). Detector banding (horizontal) was also 
noted in band 1. The Landsat 2 data evidenced banding in bands 
1, 2, and 3. The banding effects, both horizontal and diagonal, 
can be removed from all data sets using commonly available 
detector normalization and/or frequency filtering algorithms, and 
need not adversely impact use of the data for change detection. 
The banding did appear more severe in the Landsat 2 and 3 data, 
with interdetector banding appearing to be completely corrected 
for interdetector histogram equalization processing applied by 
Goddard to the Landsat 4 MSS data. We feel further conclusions 
of the radiometric and geometric properties of each sensor based 
upon sca^tergram and image differencing can not be made until the 
additional scene pairs are received and compared. 
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COMPARISON OF LAND COVER INFORMATION FROM 
LANDSAT MSS AND AIRBORNE TMS FOR HYDROLOGICAL 
APPLICATION'’ 

Janette C. Gervln 
NASA/Goddard Space Flight Center 

Yun Chi Lu 
Richard F. M*<rcell 
Computer Sciences Corporation/ 

Goddard Space Flight Center 

Greenbelt, MD 20771, U.S.A. 

ABSTRACT 

In a cooperative program with the US Amy Corps of Engineers, NASA 
is evaluating the capabilities of Landsat-4 Thematic Mapper (TM) data 
for environmental and hydrologic applications. Several sites already 
under study by the Corps were selected for this program. This 
paper will report results for one of these sites, the Clinton River 
Basin in Michigan. 

Land cover information for the Clinton River Basin derived from Landsat 
Multlspectral. Scanner (MSS) data was compared with that from airborne 
Thematic Mapper Simulator (TMS) to investigate the probable capabilities 
of the Thematic Mapper (TM) launched 'board Landsat-4 in July 1982. 

This paper reports the findings for two 7.5 minute topographic maps. 
Significant Improvements in land cover classification accuracy were 
obtained using TMS data as compared with MSS data. Overall mapping 
accuracy increased from 44 to 58 percent with an improvement from 
62 to 82 percent in the residential cetegory. A combination of four 
bands with one band in each major region of the spectrum (visible, 
near IR, middle IR and thermal IR) provided nearly as good a dis- 
crimination of land cover as all seven bands. Based on the improved 
land cove , classification accuracy of TMS, TM data has the poten- 
tial to piovide more useful and effective input to US Army Corps of 
Engineers flood forecasting and flood damage prediction/assessment 
. ode Is . 
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LANDSAT-4 MSS RADIOMETRIC CORRECTION: 
METHODS AND RESULTS 

A. Singh, PhD. 

General Electric Space Division 
Lanham Center Operations 
4701 Forbes Boulevard, Lanham, MD 2C706 


ABSTRACT 


Radiometric correction for the Landsat-4 MSS system consists of two 
parts. The first is the prelaunch calibration * the Internal lamps. 
The sensors are exposed to the known radiances or the integrating 
sphere, which allows for the calibration of the lamps. The results 
of this process are the regression coefficients (C.'s & D.'s). 

These coefficients are used in the second part - the on-line correc- 
tion of the data. In this step the detector voltage samples are 
converted into values representing the Input radiance. The step 
consists of estimating Gains and Biases of the detector transfer 
curves, modifying them using scene content information, and generat- 
ing the Radiometric Look Up Tables (RLUTs). These tables are then 
applied to the raw data (voltage samples) to yield the input radiance 
estimates (archival tape data). 


This presentation will review the on-line correction process, and 
the fine tuning of the parameters that are necessary to accomodate 
for sensor drifts. It will also contain a chronological summary of 
the parameter updates that were made during the first six months of 
operation. Finally, it will examine the Radiometric Quality Indica- 
tors (RQI) that were developed to evaluate the radiometric correction 
process, and present performance results. 
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LS-4 MSS GEOMETRIC CORRECTION: 
METHODS AND RESULTS 


J. Brooks, PhD. and E. Kimmer, PhD. 

General Electric Space Division 
Lanham Center Operations 
4701 Forbes Boulevard, Lanham, MD 2 0706 


ABSTRACT 

LS-4 MSS Geometric Correction data are generated in a two step process. 

The first step utilizes models of the spacecraft-scanner-spinning earth 
system, with smoothed ephemeris and attitude data as input, to deduce 
systematic correction data. The systematic correction data therefore 
suffer from random system pointing and spacecraft location errors. If 
control points are available for the subject scene, the second step - 
automatic correlation of control point chips to control point neighbor- 
hoods - is performed to develop input to the MSS control point location 
error filter. The output of the filter is the state error vector, which 
is used to upgrade the systematic correction data to geodetic correction 
data . 

An important offline subsystem of the LS-4 MSS Image Generation Facility 
is the Control Point Library Build (CPLB) subsystem. This hardware and 
software complex is used to generate control point chips and support 
data for use in online upqrade of correction data. Control point chips 
that represent ground truth are called geodetic or supplemental, depending 
on whether they are used to establish the model surface (with the same 
filter that Is used online) or are selected from the fully corrected 
imagery. Relative control points are designated in systematically cor- 
rected imagery, have no ground truth and can therefore be used only for 
temporal registration. 

During the last six weeks before the NASA/NOAA turnover of the system, 
an intensive effort was carried out to fine-tune both software and data 
base and to evaluate the geometric performance of the calibrated system. 
The online processing and the CPLB process provided all the tools neces- 
sary for the verification, calibration (including initial scan profile 
correction) and evaluation of LS-4 geometric processing. In this presen- 
tation, we will describe the methods and results of that effort, including 
changes to system and data base that will affect the LS-4/EDC products. 
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Spectral Characterization of the Landsat-4 Thematic Mapper Sensors 

by 

Brian L. Markham and John L. Barker 
NASA/Goddard Space Plight Center 
Earth Resources Branch/Code 923 
Greenbelt, MD 20T71 


A summary and analysis of data collected by Hughes/SBRC* on the spectral characteristics 
of the Landsat-4 and Landsat-4 backup Thematic Mapper instruments, the protoflight 
(PF) and flight (F) models, respectively, is presented. Tests were conducted 
on the instruments and their components to determine compliance with two sets 
of spectral specifications: Band-by-band spectral coverage and channel-by- 

channel within-band spectral matching. 

Spectral coverage specifications were placed on: (1) band edges — points at 

50% of peak response, (2) band edge slopes — steepness of rise and fall-off of 
response, (3) spectral flatness — evenness of response between edges, and (4) 
spurious response — ratio of out-of-band response to in-band response. Compliance 
with the spectral coverage specifications was determined by analysis of spectral 
measurements on the individual components contributing to the spectral response: 
filters, detectors, and optical surfaces. Th'. protoflight and flight model 
TM * s used filter pieces cut from the same substrate and detectors from the 
same batch (except hand 6) ; any differences between the calculated relative 
spectral responses (RSR) resulted from optics differences (except band 6). 

The RSR 1 s for the reflective bands were similar between PF and F except for 
the within-band flatnesses. Both the PF and F reflective band calculated 
spectral responses were within specifications except for bands 2, 3 and 7 
spectral flatnesses and the band 5 upper band edge: specification, 1750 +20nm; 
calculated as 1784 nm. In band 6, the PF and F showed fundamentally different 
spectral responses. The PF upper-band edge was detector determined at 11.66 pm 
(varying with temperature) ; the F upper band edge was filtered determined at 
12.43 pm. The PF band 6 was spectrally out-of-specification, but its detectors' 
sensitivities were better than th 2 F, which was basically within specifications 
spectrally. 

Out-of-band responses for all bands were within specification. Bands 1 and 3 
had some sensitivity to near-IR radiation. Band 1 filters had transmission 
peaks at 800 and 885 nm of 0.5% and 0.7% respectively. When measured on the 
flight model, an approximately 1 count contribution to band 1 resulted when 
the radiance between 776 and 905 nm resulted in 100 counts in band 4. Band 3 
filters had peaks in out-of-band transmission at 945 and 1000 nm, of 2.8% and 
1.2% respectively. The impact of this on the band 3 response has not been 
determined . 


'Data and assistance in interpretation were provided by Hughes/SBRC personnel: 
Richard W. Cline — spectral coverage determination 
Jack C. Lansing — spectral matching 

Donald G. Brandshaft — out-of-band response, light-leaks 


The spectral matching specification stated that "after system calibration, the 
peak-to-peak signal variations between channels within any of the first five 
bands and band seven, when all channels are viewing the same scene radiance, 
shall be less than 0.5 percent of the minimum saturation levels for the two 
test conditions whose parameters are given in . . .(a linearly varying spectral 
radiance and a flat spectral radiance)." Initial test plans called for using 
channel-by-channel relative spectral response curves to determine compliance 
with the spectral matching specification. As channel-by-channel relative 
spectral responses were not measured for the TM's, an alternative test was 
devised to assess the spectral mismatch between channel in a band. This test 
invoked calibrating the individual channels on the 1.2 meter integrating 
sphere and then recording the mismatch in their outputs to a spectrally different 
source, the TM calibrator (modified by filters). The protoflight teat gave 
out-of-specification results which appeared to be attributable to spatial non- 
uniformity of the calibrator source. A refined test was used for the flight 
model testing, using the calibrator source with and without filters for the 
two targets. With the exception of band 4, which showed a 1.7% mismatch, all 
bands were within specifications. Calculations using the relative spec .’al 
response data for the 5 MSS sensors 'MSS 1, 2, 3, 4 and 4-backup) showed that 
the flight model TM had comparable or better spectral matching than the MSS 
sensors. 

An examination of the along-scan line spread function for the flight model TM 
revealed several minor white light leaks in the primary focal plane bands (1- 
4)~ i.e., light getting to the detectors without passing through the spectral 
filters. The magnitude of these light leaks is dependent on the spectral 
character of the illumination. For the odd channels of band 1 (magnitude of 
the light leaks comparable for all detectors in a half-band), which was the 
worst observed, with the TM calibrator 'white' light source, a light leak at 
13.1 IFOV off the detector center made about a 1» contribution to the signal. 

The location and shape of the light leaks suggests they are associated with 
the slots at the sides of the individual band assemblies. It is believed the 
protoflight model has comparable light leaks. 
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Prelaunch Absoluce Radiometric Calibration of the Reflective Banda on 
Landsat-4 Protoflight Thematic Mapper 


J. L. Barker (NASA/GSFC) , D. L. Ball and K. C. Leung (CSC) 
and J. A. Walker (SBRC) 


Abstract 

Landsat-4 Early Results Symposium '23 Feb 83) 


We present results on the absolute radiometric calibration of the 
Landsat-4 Thematic Mapper (TM) as determined during prelaunch tests 
with a 122-cm Integrating Sphere (IS). Discussions of the radiometric 
requirements specified for the TM reflective bands, the equipment 
used to calibrate the instrument and the procedures used are included. 
Unlike the photomultiplier tubes used on the Multispectral Scanner, 
the solid state detectors on the TM an linear in their response to 
radiance. The sample response shown in the accompanying figure is 
linear to 0.4 percent of full scale in the standard error of the 
estimate for a linear fit. 

In the table below we summarize results for the dynamic range, 
average gain and offset and signal-to-noise ratio (SNR) for the TM 


reflective bands 



ORIGINAL PAGE IS 
OF POOR QUALITY 



Dynamic 

Range 

Band 

RMIN 

RMAX 


(mWcm 

-1 -1. 
sr pm ) 

1 

-0.152 

15.842 

2 

-0.284 

30.817 

3 

-0.117 

23.463 

4 

-0.151 

22.432 

5 

-0.037 

3.242 

7 

-0.015 

1.700 


Gain 

- 2-1 

(counts/ (mWcm sr pm 

15. 7810.543!: 
8.1010.892 
10.6?l0. 76% 
10.90H .232 
77.2410.532 
147.1210.722 


Signal-to-Noise 

Ratio 


Of f set 

Spec 

Obs 

(counts) 

(MSL) 3 

(MSL) 3 

2.5810.18 

85 

152192 

2.4410.16 

170 

2811172 

1.5810. 20 

143 

235162 

1.9U0.22 

240 

341192 

3.0210.11 

75 

180162 

2.4U0. 20 

45 

1751142 


1RMIN is the largest minimum radiance observed in the band (0 counts). 
^RMAX is the smallest maximum radiance observed in the band (255 counts). 
*MSL is the minimum saturation level radiance for the band. 


Three detectors failed to pass the minimum SNR specified for their respective 
bands; band 5, channel 3 (dead), band 2 channels 2 and 4 (noisy or slow re- 
sponse) . 

Estimates of the absolute calibration accuracy for the TM show that the 
detectors are typically calibrated to five percent absolute error for the 
reflective bands; ten percent full scale accuracy was specified. Ten tests 
performed to transfer the detector absolute calibration to the Internal 
calibrator show a five percent range ar full scale in the transfer cali- 
bration; however band 5 in two cases showed a ten percent difference and a 
seven percent difference. These errors are given under the assumption of 
identical operating conditions for all transfer tests; known temperature 
sensitivities fot .e TM have not yet been included. 
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Characterization of Radiometric Calibration of LANDSAT-4 

TM Reflective Bands 

J. L. Barker, R. B. Abrams, D. Ball, and K. C. Leung 


Knowledge of the absolute radiometric accuracy of Thematic 
Mapper (TM) LANDSAT-4 data is critical for arithmetic 
spectral transforms such as those used to determine atmos- 
pheric correction path radiance. In this paper the radio- 
metric performance of the Landsat-4 TM is characterized. 
Radiometric calibration using nominal gains and offsets 
versus using the internal calibration system (IC) is dis- 
cussed. Prelaunch and postlaunch IC, image and background 
data are examined. Early results indicate that IC lamp 2 has 
changed relative to prelaunch. For any given band, the change 
in counts relative to November 2, 1982, for each single-lamp 
IC state is the same for all three states within .5 percent. 
When prelaunch March 9, 19>2 data is used as a reference, 
lamps 1 and 3 have the same change in counts wxthin .8 percent 
while lamp 2 diflers from lamp 1, on the average, by 2.2 
percent. All but two channels have performed normally. 

Band 5, channel 3 and band 2, channel 4 have been replaced 
with data from band 5, channel 4, and band 2, channel 5 
respectively. Gain changes relative to November 2, 1982, 
for channels within a band vary as a group within .5 percent. 
Change in gain of primary focal plane channels was 1 percent 
for bands 2 and 4 and 2 percent for bands 1 and 2 from 
August to October, and less than .8 percent from October 
through December (see Table 1) . The variability in the cold 
focal plane of up to 7 percent in band 5 and 4.5 percent in 
band 7 indicates a need for IC correction. The gain in 
band 1/channel 1 was under valued by approximately 5 percent 
on all data processed prior to December 22, 1982, due to 
clipping of the calibration pulse in the ground-processing 
system. 

Background noise ranges from .5 to 1.27 digital counts 
depending on the band (see Table 2) . Noise in image data 
from flooding lamp (prelaunch calibration) and water scenes 
is slightly higher, .5 to 1.7 counts 'see Table 3). Discrete 
fourier transforms of data from all 96 reflective channels 
(first scan) in a scene from Memphis, Tennessee, August 22, 
1982, show coherent noise at 32.7 kHz in band 1/channels 2, 

3, 4, 5, 6, 8, 14, and 16; band 3/channels 4 and 8, and 
band 4/channels 4 and 5 with peak-to-peak amplitudes 
ranging from .7 to 1.4. 



Possible improvements in calibration processing are pre- 
sented. Forward and reverse average differences for back- 
ground flooding lamp and water scenes are erratic up to one 
digital count. This might indicate a need for recalibration 
of forward and reverse prelaunch absolute radiometry. Preci- 
sion is improved typically by 50 percent by increasing the 
pulse integration width from 31 to 31-41 minor frames 
depending on the band. Dropping 4 minor frames preceeding 
a transition improves precision for cold focal plane bands 
by 35 percent. 


Table 1 

Postlaunch Change in Gain in Parts Per Thousand Relative 
to March 9, 1982 Prelaunch Values 
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THERMAL BAND CHARACTERIZATION OF LANDSAT-4 
THEMATIC MAPPER 

Jack C. Lansing/Santa Barbara Research Center 

and 

John L. Barker/Goddard Space Flight Center 
ABSTRACT 

Observations With Time 

A quick-look monitor In the spacecraft control center wil used to 
measure the TM Band 6 shutter background and the 34.7°C Internal block 
body signal on over 50 dates. Variability of the shutter background 
temperature varied from 7 to 11°C. For ten specific images the digital 
counts of the calibration data were measured. The average pulse value 
of the block body peak, CB, decreased from 174 to 149 counts while 
the shutter background counts, CS, varied as a function of shutter 
temperature, TS, from 77 to 85 (CS * 3.05 TS(°C) + 53.27). The net 
calibration peak, CB-CS, decreased from 95 to 64 between August 22nd 
and December 20th. Relative internal gains between the four channels 
were calculated and compared to pre-launch values; they showed changes 
over 9 months of up to 5%. 512 x 512 subsections of the original ten 

daytime scenes showed target counts, TS, that ranged from 135 down co 
62, with a range of standard deviations from as low as ± 0.5 for a pure 
water scene off Bostor to as high as ± 4.1 for an August 22nd scene 
over Arkansas. 

Image Destripinq 

Frequency histograms of numbers of pixels vs. digital counts from 
a night scene of the Buffalo area of 22 August 1982 were used to deter- 
mine channel gain relative to the mean and to discern a systematic 
along-scan pattern in a difference between forward and reverse scan 
counts of up to 0.5. These results were used to produce a corrected 
digital image, which, in turn, exposed minor adjustments necessary to 
further correct the channel gains and offsets. Individual gains and 
offsets were calculated for the four channels. The final values may 
be used to produce a destriped image of Band 6. 



At-Satellite Radiance 


An Illustrative thermal band radiometric calibration graph was 
proiuced to relate knjwn spectral radiance to digital counts. Calibration 
lines were calculated for Initial turn on of the cold focal plane In August 
through the 60% gain loss in Band 6 over 5 months until outgassing of the 
cooler In early January, when the gain was restored to its original 
value. 

The specified range of 260K to 320K corresponds to a range of 
approximately 50 to 195 counts, which is centered In the 0 to 255 
counts available. 

Noise Equivalent Temperature 

An area of 600 120m pixels in Lake Erie was used to calculate noise 
equivalent temperature difference, which was 0.10 K at 300 K, neglecting 
the forward to reverse scan difference. 

Target Temperature After Atmospheric Correction 

The calibration data and the Buffalo scene with the corrections 
mentioned and with preliminary estimates of the atmospheric transmission 
and radiance were used to make a temperature estimate for an area of 
Lake Erie, of 18.5 to 19.9°C. Local records of the temperature showed 
21°C. 


Landsat-4 TM Thermal Band 6 Radiometry 

Channel 4 

260 K (-13°C) 320 K (47°C) 
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PRELIMINARY EVALUATION OF THE 
RADIOMETRIC CALIBRATION OF LANDSAT 4 THEMATIC MAPPER DATA 
BY THE CANADA CENTRE FOR REMOTE SENSING 


AUTHORS: J- MURPHY W. PARK. A. FITZGERALD 


PRESENTED BY JENNIFER MURPHY 


CANADA CENTRE FOR REMOTE SENSING 
2464 SHEFFIELD ROAD 
OTTAWA, ONTARIO 
CANADA 


ABSTRACT 

The technique being evaluated by the Canada Centre for 
Remote Sensing (CCRS) fpr the radiometric correction of 
Landsat 4 Thematic Mapper ( TM ) data is discussed. 
Preliminary resuits on the removal of radiometric striping, 
caused by inequailties in the calibration c r individual 
detectors within each band, are presented. The destriping 
method was originally developed by CCRS for processing 
Landsat Mu 1 1 1 s pe c t r a 1 Scanner (MSS) data and uses the scene 
histograms to equalize the responses oi the individual 
detectors within each band. For the TM case, it may be 
applied to the lorward and reverse scans either separately 
or in combination CCRS is also evaluating the radiometric 
calibration of test scenes using the absolute pre-launch 
calibration tables supplied by NASA. Quantitative estimates 
of residual striping in the r ad i ome t r ic a 1 1 y corrected TM 
data are provided- 
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Evaluation of the Radiometric Quality 
of the TM Data Using Clustering and 
Multlspectral Distance Measures 

L.A. Bartolucci, M.E. Dean and P.E. Anuta 
Laboratory for Applications of Remote Sensing 
Purdue University 
West Lafayette, Indiana 


Int ro d uct Ion 

The primary objective of this investigation was to evaluate the 
radiometric quality of the Thematic Mapper (TM) data for classification 
and identification of earth surface features. In addition, a comparison 
of the information content between TM and MSS data sets was carried out. 
Finally, the TM data were utilized to map the thermal effluent discharge 
into a river ecosystem from a nuclear power plant, an application only 
possible until now through the acquisition of thermal infrared scanner 
uata from aircraft rltitades. 


Procedures 

Padiomet rical ly corrected (A-tape) and geometrically corrected (P-tape) 
TM data from three different geographic locations were utilized in this 
investigation. A description of these data sets is given in Table 1. 


Table 1 . TM Data Seta 


Scene ID 
40037-16031 
40049-16264 
40101-16025 


Date 

8 / 22/82 

9 / 03/82 

10 / 25/82 


Location 

Arkansas 

Iowa 

Illinois 


The Initial examination of the TM data consisted of an Inspection of 
the histograms for each band to determine the dynamic range of the data, 
the shape of the distributions, arid to verify whether empty bins were 
Introduced by the radiometric correction process. 

In order to assess the effect of the geometric correction process 
(resampling effects) on the radlometry of the new resampled pixels, the 
means and variances for homogeneous and heterogeneous areas on th«t ground 
scene were calculated using both the A and P tape data sets. Furthermore, 
to determine the Influence of resampling on the structure of the data, 
clustering of the A and P tape data sets of the same areas on the ground 
were performed. The resulting cluster classes from the original and 
resampled data sets were then compared using a transformed divergence 
measure. The results of this multispectral distance measure were also 
used as a feature selection criterion to obtain the optimum combinations 
of spectral bands for dimensionality reduction purposes. 

The intrinsic dimensionality of the TM data was also studied by 
means of a principal components analysis. A Karhunen-Loeve transformation 
va 3 applied to both the MSS A and P and to the TM P data tapes of the 
Chicago O'Hare airport test site. Statistics used to generate the 
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transformation matrices for each of these date sets were calculated from 
samples of the data taken from every 5th line and 5th column In the scene. 

The resulting transformed components were scaled to be proportional to the 
square root of their respective eigenvalues, with the first component 
set equal to a range of 0*255. 

In multidimensional space, such transformations will not alter the 
Euclidean distance between data points. However, If the resulting 
components are given some relative magnitude or length other than values 
which are proportional to the square roots of their respective eigenvalues, 
then the Euclidean distance between data points will be altered. Clustering 
algorithms are often highly sensitive to such movements or changes In 
position of data points with respect to each other. Therefore, It Is 
Important for clustering and classification purposes that the resulting 
transformed components of such orthogonal transformations be (each) correctly 
scaled . 

To assess the utility of high (spatial) resolution thermal infrared 
data for thermal pollution studies., tm .*at>. over an area south of Chicago, 
Illinois, were classified using a layered clas3lfler to first discriminate 
water from everything else with the aid of the reflective infrared binds, 
and then to perform a level slicing of the thermal band for only che 
resulting water pixels, thus producing a temperature map of the thermal 
plume. The linearity of the calibration of the thermal bar J was also 
studied for various temperature ranges. 
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Results 

An evaluation of the histograms from the A and P tape TM data sets 
shows that the data in the three reflective Infrared bands have a larger 
dynamic range than the data in the visible and thermal bands. It was also 
observed that the shape of tha distributions in the three visible bands, 
especially in the blue band, deviated considerably from a normal curve, 
approximating a distribution highly skewed towards the higher response values. 
This phenomenon was not observed in the Landsat 4 MSS bands corresponding 
to the visible portion of the spectrum. This histogram also showed that 
the A tape data from both the MSS and TM systems contained a number of 
empty bins. As it would be expected, empty bins were not found in either 
the MSS or TM P tapes, since the gray level interpolation applied to the 
geometrically corrected data would eliminate such effects. 

From past experience at LARS, there was concern about the magnitude of 
the radiometric changes that would occur on resampled (gray level interpolated) 
pixels. To determine quantitatively the effect of the interpolation, 
individual pixels and sets of pixels representing the same features on the 
ground were examined. The results indicate that the interpolation did not 
affect the means and variances of the spectral response of ground cover types. 
In fact, the smoothing effect that the cubic convolution was expected to 
have on the radiometry of the data was not evident in these results as 
Illustrated in Table 2. 

In order to further evaluate the effects of the cubic convolution on 
the structure of the geometrically corrected TM data, the same areas on the 
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ground were selected In the A and P TM Images and both data sets were 
clustered. The statistical parameters (means and covariance matrices) that 
describe the resulting cluster classes from both the A and P TM data sets 
were found to be virtually Identical; but to determine more quantitatively 
the similarity between corresponding class pairs, a Transformed Divergence 
(Spectral Separability) algorithm was used. Table 3 Illustrates the results 
of this separability analysis. The Transformed Divergence values range 
from zero (Identical spectral classes) to 2000 (completely different spectral 
classes) . 

The same spectral separability measure was also used as a feature 
selection function In order to determine the best combination (subset) of 
bands that would provide the most accurate multlspectral classification 
results. Also, the MSS data for the same areas on the scene were clustered 
and the results were compared with the resulting cluster classes from the 
TM data In an attempt to learn more about the Information content from the 
two multlspectral scanner systems. 

Finally, a temperature map of a nuclear power plant cooling pond and 
the thermal plume discharged Into the Illinois river was produced using a 
Layered Classifier and calibrating the relative digital counts in the 
thermal IR band Into degrees Celclus. It will also be shown that the 
relationship between the thermal IR radiation and the corresponding radiant 
temperature Is not a linear function for the 10.4-12.5 vim band and for the 
range of temperatures from 260K to 320K. 


Classes Considered Class Pairs D (Combination of all bands) 
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ANALYSIS AND EVALUATION OF THE LANDSAT-4 MSS AND TM SENSORS AND 
GROUND DATA PROCESSING SYSTEMS-EARLY RESULTS 

Ralph Bernstein 
Jeffrey B. Lotsplech 
IBM Palo Alto Scientific Center 
1530 Page Mill Road 
Palo Alto, CA 94303 0821 USA 
(415) 855-3126 


ABSTRACT 

This Investigation deals with the analysis of the new Landsat-4 sensor 
data and ground data processing. The analysis involves the assessment of 
the data radiometric calibration accuracy, geometric correction accuracy, 
data resolution and data quality. In addition. Image research dealing with 
advanced Information extraction techniques are planned. Both Landsat-4 
Multlspectral Scanner and Thematic Mapper data have been acquired and 
processed. The data processed to date Include both raw and processed data 
from three geographic areas of Interest. Several digital Image processing 
techniques have been applied to the data to Improve data radlcmetry and 
quality. 

In the area of radlometry, although the quality of the data Is very good. 
Incomplete radiometric calibration processing Is evident from the data 
statistics and image products. Improved radiometric calibration results 
have been developed, and consistently demonstrate less than one quantum 
level of relative radiometric error for all detectors signals. This 
processing has substantially reduced the striping evident In the original 
data. To assess the Information content and quality of the data, entropy 
analysis and principal component processing of the data have been performed. 
Also Intensity, hue, and saturation data presentation experiments have 
been conducted. The data exhibit excellent information potential, although 
data entropy studies show that the full range of the sensor performance has 
not yet been achieved. In the area of data quality, some problems are 
evident, due to failed detector compensation. Improved techniques have 
been developed to compensate for these failures, and the results are 
presented. These show that better techniques than the line repeat 
procedure should be used. Some sub-pixel band-to-band misregistration is 
evident, and registration statistics are presented. 

Resolution studies using subimage areas of well defined features were 
performed using precise ground truth data. In addition, data edge analyses 
studies have been conducted. These results are presented In image and 
statistical form, ar.d have shown that the TM data exhibits very good data 
quality and reliability, and can be used for accurate mensuratioi studies 
at a suitable scale. 
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In -progress Absolute Radiometric 
Inflight Cslibration of the Landsat-4 Sensors 

K. Castle (Universi* j of Arizona), M. Dinguirard (CERT), C. E. Ezra (USDA), R. G, 
Holm (UA), K. J. Jackson (USDA), C. J. Kastner (UA), J. M. Palmer (DA), R. Savage 
(ASL), and P. N. Slater (DA). 

Abstract 

We describe our approach to providing periodic inflight absolute 
radiometric calibrations of the Landsat-4 sensors by reference to selected, 
instrumented ground areas, and we present the results of some early ground 
measurements and computer simulations. 

In pursuing such a calibration approach, the first step is the selection 
of a suitable ground reference site. The requirements are that the site should 
contain areas that are ( 1 ) large, of the order of several square kilometers to 
minimize atmospherically induced adjacency effects; ( 2 ) flat, diffuse and of 
uniform relectance, to minimize the difficulty of accurately mapping surface 
spectral radiance variations; (3) of different reflectances, to provide several 
points on the calibration curve; (4) in a region where the sky is often clear 
and the aerosol content of the atmosphere is low; (5) within convenient 
travelling distance for the measurement team. The White Sands Missile Range 
in New Mexico is such a site, containing an extensive white "alkali flats” area 
of gypsum and other areas of dark soil and water. 

The second step is the accurate measurement of the spectral reflectance 
of the selected ground areas and the determination of certain atmospheric 
characteristics (the Rayleigh, Mie, O 3 spectral extinction coefficients, and h 2 & 
vapor concentration) during the morning of the sensor overpass. This involves 
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Che accurate absolute calibration ot the instruments used and the aerial 
mapping ot spectral variations over the selected ground areas. In preliminary 
measurements we are using a manually operated solar radiometer, two 
transpot. table Barnes Model 1200 Mhk radiometers, an Exotech Model 10CA, and an 
RMk-10 Radiometer. 

The third step is the reduction of the measured data and the use of these 
data in an appropriate atmospheric radiative transfer program to predict the 
radiance levels at the entrance pupil of the sensors across each of their 
spectral bands. 

The last step is the comparison of the radiance level data with the 
digital counts for the images of the selected ground areas. This comparison 
will be made first for the raw data in CCT-D form, then for data to which the 
internal calibration and the histogram equalization procedure has been applied 
(CCT-A) and finally to the radiometrically corrected and cubic convolution 
resampled data in CCT--P form. 

At this writing, attempts have beer made to collect data at White Sands 
on 18 December 1982, 3 and 19 January 1983, and A February 1983. Because of 
cloud cover, only data for 3 January 1983 were collected. On that date the 
ground had a 80 mm snow cover which provided a very uniform high reflectance 
surface in bands 1 through 4 and a low uniform reflectance in bands 5 and 7. 
These data have been reduced together with field and laboratory measurements 
of gypsum and are presented at this symposium together with Langley plot data. 

These preliminary measurements are being made partly as an aid to 
defining the characteristics of field equipment to be constructed ana 
calibrated for later use on this program. Computer simulations are also being 
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made to predict the periods of the year for which the various TM bauds will be 
unsaturated by the radiance of the main gypsum reference area at White Sands 
and to determine the accuracy to which the various ground-based measurements 
have to be made. Some of these results are presented at this symposium. 
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Landsat-4 Thematic Mapper Calibration and Atmospheric Correction 

Warren A. Hovls 

National Oceanic and Atmospheric Administration 
National Environmental Satellite, Data, and Information Service 

Washington, D. C. 20233 

ABSTRACT 

Ti'e Landsat-4 Thematic Mapper, with Its wide spectral coverage 
and digitization to 8 bits per word. Is a large step forward In the 
direction of quantitative radlometry from the Mu’tl Spectral Scanner 
(MSS) flown on all four Landsats launched to date. In order to uti- 
lize the quantitative accuracy built Into the Thematic Mapper effec- 
tively, more attention must be paid to calibration before launch, 
changes of calibration with time In orbit, and atmospheric Interfer- 
ence with the measurements, especially In the 450 to 520 nanometer 
band. All of these factors are Important if we are to realize the 
goal of such an Instrument; namely, determining the upwelled radiance 
that would be measured at the surface. Recent experience with the 
Coastal Zone Color Scanner (CZCS) program has led to procedures where- 
in Rayle'-jh correction factors can be generated utilizing simultaneous 
surface truth data that empirically give correct upwelled surface 
radiances, despite errors In sensor calibration, solar spectral Irradl- 
ance measurements, and reported values of Rayleigh optical depth. These 
techniques also offer sensitive tests for change in calibration, 
especially at shorter wavelengths. Previously launched Instruments, 
such as the CZCS, have shown that calibration changes first and to the 
largest degree at the shorter wavelengths, with lesser changes as wave- 
length Increases. These techniques, applied to Thematic Mapper data, 
are utilized to calculate a Rayleigh correction factor that, together 
with geometric terms, will give an accurate correction for this portion 
of the atmospheric contribution to the signal. Long term observation . 
of the Rayleigh corrected radiance over clear water will be a sensitive 
Indicator of any change In calibration as the sensor ages In orbit. 
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SCAN-ANGLE AND DETECTOR EFFECTS IN 
THEMATIC-MAPPER RADIOMETRY 


Michael D. Metzler and William A. Ma 111a 

Environmental Research Institute of Michigan 
Ann Arbor, Michigan 48107 


S unwary 


Objective 

The objective of ojr investigation is to quantify the performance 
of the Thematic Maoper (TM) as manifested by the quality of its image 
data, in order to suqqest possible improvements for data production and 
assess the effects of data quality on its utility for land resources 
applications. The major emphasis of our work thusfar has been on the 
radiometric characteristics of TM data. 


Approach 

Computer-compatible tapes of raw data (CCT-B), radiometrically 
corrected data (CCT-A), and geometrically corrected data (CCT-P) have 
been examined for two scenes, 40049-16262 (Iowa) and 40037-16031 
(Arkansas), with concentration on the CCT-B and CCT-A data. Histo- 
grams, scan-line (across-track) and along-track averages, Fourier 
transformations, and gray maps have been produced and analyzed, for 
each band, each detector, and each scan direction. 


Initial Findings 

Our initial examination of Thematic Mapper data yielded an overall 
positive reaction to its quality and information content. We did, how- 
ever, observe and quantify some relatively low-level artifacts relatable 
to the sensor and data processing stages. We also have noted and quanti- 
fied scene-related scan-angle effects that could be important in analysis 
of full scenes and in comparing data from different portions of a scene. 

Horizontal stripes and discernable banding are common artifacts in 
image data from sensors that scan arrays of detectors. We found evidence 
of striping in corrected TM data due to residual between-detector cali- 
bration differences and quantization effects like those found in MSS data. 
Figure 1 illustrates one pattern of differences found in Band 2 between 
averages of successive scan lines across a homogeneous scene. Effects 
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of unequal analog-to-dlgital conversion bln sizes in CCT-6 data and empty 
quantization levels in CCT-A data are largely masked by the interpolation 
process used to produce CCT-P data. 

In addition, we detected a new type of banding caused by the com- 
bination of bidirectional scanning by the 16-detector arrays and a sys- 
tematic droop of signal values during the active portion of each scan. 

The resulting differences between forward and reverse scans are Illus- 
trated in Flqure 2, as a function of scan angle for TM Band 1. This 
effect produces swaths or bands of differing signal levels that are most 
pronounced near the frame edges and minimal at the frame center. Similar 
behavior was detected In Bands 2, 3, and 4, relatively little effect In 
Bands 5 and 7, and a substantially different behavior <n the thermal band 
(Band 6). The markedly different full-frame averages for forward and 
reverse scans In the thermal band are shown In Figure 3. 

Tne scan -related banding In the reflective bands was superimposed 
on substantial scene -related scan-angle effects In the data. Again for 
Band 1 , the full-frame averages for forward and reverse scans are pre- 
sented In Figure 4. The dominant scan-angle effect here Is due to dif- 
ferential backscatterlng by the atmosphere, differing path lengths, and 
bidirectional reflectance effects on the ground. It was expected to be 
greater than seen before In Landsat MSS data due to the larger scan 
angles for all bands and the shorter wavelength of Band 1. Differences 
were observed between bands due to spectral phenomena. 

An empirical first-order model to correct for the scan-direction 
banding effect was derived for the reflective bands. Refinement of the 
model awaits further understanding of the sensor phenomena driving the 
effect and other scene data for validation. We also recommend that 
procedures be investigated to normalize the scene-related scan-angle 
effects, as has been done previously for aircraft and Landsat MSS data. 
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'IGURE 1. DOWN-TRACK MEAN SIGNAL, BAND 2, 
ILLUSTRATING DETECTOR 
CALIBRATION DIFFERENCES 



FIGURE 2. DIFFERENCE IN CROSS-TRACK 
MEAN SIGNAL BETWEEN FORWARD 
AND REVERSE SCANS, BAND 1 


CCT-AT 4-0049-16-262 TVI P027 ROJO 
BANC 6 DETECTOR 4 



CCT-AT 4-0040-16-262 TM P027 ROJO 



FIGURE 3. CROSS-TRACK MEAN SIGNAL, BAND 6, 
ILLUSTRATING SCAN DIRECTION 
EFFECT 


FIGURE 4. ILLUSTRATION 0^ SCAN-ANGLE 
EFFECT, ACROSS A FRAME 
(BAND 1) 
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THEMATIC MAPPER SPECTRAL DIMENSIONALITY AND DATA STRUCTURE* 

E. P. Crist and R. C. Clcone 

Environmental Research Institute of Michigan 
Ann Arbor, Michigan 48107 


Abstract 


Thematic Mapper data, simulated from field and laboratory 
spectrometer measure».'.ents of a variety of agricultural crops and a 
wide ronge of soils, are analyzed to determine their dispersion In the 
six-space defined by the reflective TM bands (l.e., excluding the 
thermal band). While similar analyses of MSS data from agricultural 
scenes have found that the vast majority of the MSS data occupy a 
single plane, the simulated TM data are shown to primarily occupy 
three dimensions, defining two Intersecting planes and a zone of 
transition between the two (see Figure). Viewing the "Plane of 
Vegetatlor” head-on provides a projection equivalent to the single 
plane of MSS data. The "Plane of Soils" and transition zone represent 
new Information made available largely as a result of the 
middle- Infrared bands Included on the Thematic Mapper. Variation In 
this transformed data space Is described In terms of the physical 
properties of the vegetation and soils In the data set, and the 
potential utility of the new Information sources Is highlighted. 


*The work described herein was supported by the Inventory Technology 
Development project of AgRISTARS under NASA contract NAS9-16538. 
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MTF Analysis of Landsat-4 Thematic Mapper 
Robert Schowengerdt 
University of Arizona 


A research program to measure the Landsat-4 Thematic Mapper (TM) modula- 
tion transfer function (MTF) is described. The research is being conducted by 
the University of Arizona arid NASA Ames Research Center under the contract "An 
Investigation of Several Aspects of Landsat-D Data Quality." 

Measurement of a satellite sensor's MTF requires the use of a calibrated 
ground target, i.e. the spatial radiance distribution of the target must be 
known to a resolution at least four to five times greater than that of the 
system under test. Any radiance structure smaller than this will have a small 
effect on the calculated MTF. Calibration of the target requires either the 
use of man-made special purpose targets with known properties, e.g. a small 
reflective mirror or a dark-light linear pattern such as line or edge, or use 
of relatively high resolution underflight Imagery to calibrate an arbitrary 
ground scene. Both approaches will be used in this program; in addition a 
technique that utilizes an analytical model for the scene spatial frequency 
power spectrum will be Investigated as an alternative to calibration of the 
scene. The test sites and analysis techniques to be used In this program are 
described in this paper. 
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Thematic Mapper Intraband Radiometric Performance 
Hugh H. Kleffer, Eric Ellason, Pat Chavez, 

Raymond Batson & Warren Borgeson 
U. S. Geologic Survey, Flagstaff, Arizona 86001 

In this early report, we have concentrated on a radlometrlcally 
corrected, geometrically raw (A data) and a fully processed (P data) Images of 
the Washington D.C. area (ID 40109-15140). In A-data, every 16th f iage line 
Is from a single detector. To allow recognition of features In A data, a 
first order geometric correction consisting of a 46 sample shift between 
forward and reverse scan directions was applied. The radiometric correction 
Incorporated In A-data consists of applying a Radlometlc look-up table (RLUT) 
constructed to equalize throughout a scene the mean and standard deviation of 
each detector, with the constraint of never compressing two DN levels Into 
one. The same RLUT Is used for the forward and reverse scan directions. 

For our statistical study, we used a set of lines (512 A, 1024 P) across 
the center of the scene. A sub-scene centered over Chesapeake Bay and 
approximating a flat field was used for characterizing subtle radiometric 
differences and noise problems. 

There are small differences between the average DN for the 16 detectors 
In each band, typically 0.8 DN; the standard deviations also differ typically 
by 0.4 DN. Differences of mode in a band are as large as 6 DN. Some DN 
levels appear to be strongly favored over adjacent levels. Differences 
between forward and reverse scans are approximately 0.1 DN In the mean, and 0 
to 0.6 DN in the standard deviation. Geometrically resampled images (P data). 
In which the fixed correspondance between lines and single detectors is lost, 
are statistically slmiliar to the ensentole of detectors in each band In A 
data. 

The magnitude of inter-detector variation is readily seen by making an 
image of the first derivative in the vertical (line) direction of a flat field 
and stretching progressively wider ranges to gray; most detectors differ from 
their neighbors by 1 to 2 DN. 

The effective resolution in radiance is degraded by a tendency for the TM 
to avoid certain DN levels by about a factor of two. These levels are 
consistent over all bands and detectors, and are spaced by an average of 4 
DN. In band 6, level 127 is avoided by a factor of 30. This behavior is 
masked by resampling in the P data. 

At high contrast boundaries, some of the detectors in Band 5 commonly 
over- or under-shoot by several DN and require on the order of 50 samples to 
recover; this behaviours occurs erratically. 

A coherent sinusoidal noise pattern is evident in Detector 1 of Band 3. 
One-dimensional Fourier transforms show that this "stiching" pattern has a 
period of 13.8 samples with a peak-to-peak amplitude ranging from 2 to 5 DN. 
Oscillations of the same frequency, but about half this amplitude, occur for 
one other detector In this Band. Noise with a period of 3.24 samples Is 
pronounced for most detectors in Band 1 and for 2 detectors in band 3, 
insignificant in Band 2, weak in Band 4, and obscured by noise In Bands 5 and 
7. 

A set of adjacent detectors in Band 1 has a change of response with a 
period of several scans. These detectors will have a reponse similiar to 
their neighbors for several scans, then for several scans will have higher or 
lower DN numbers than their neighbors. There is no apparent pattern to this 
drift. 

The random noise level of each detector was characterized by the standard 
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deflation of the first derivative In the sanple direction across a flat field. 
By this measure, the noise level Is: Band 1, 1.8 to 2.8 DN; Band 2, 0.74 ON 
representative, one detector each at 1.1 and 1.7; Band 3, 0.8 representative, 

4 detectors at 1.0 to 1.3; Band 4, 1.0 representative. Band 5, 1.5 to 2.1; 

Band 7, 1.8 representative, one detector at 3.1. 

The resampling algorithm used In generation of P data Incorporates the 0 
DN values which occupy non-data areas along the left and right edges of a 
frame. This causes 'fuzzy' frame edges and adds anomalous low and high DN 
values In the resulting frame histogram. In order to avoid adverse effects In 
appllctlons Involving scene statistics, such as clustering techniques, these 
values should be reset to zero after the resampling. 

A principal component analysis Indicates that a composite of the first 
component of Bands 1, 2, and 3; of Band 5 and 7; and Band 6 contain 97% of the 
Information and has reduced the effect of noise. 

A periodic noise removal algorithm developed by Chavez, et. al. has been 
adapted and applied to the coherent noise of Band 3, detector 1; the technique 
appears to woilc quite well. 

The geometric fidelity of the GSFC fllmwrlter used for Thematic Mapper 
Images was assesed by measurement with accuracy better than 3 micrometers of a 
test grid . The film output has scale errors In both the sample and line 
directions corresponding to 4 pixels In a TM Image, and skew corresponding to 
?■ pixels. 
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An evaluation of Landsat-4 Thematic Mapper data for 
their geometric and radiometric accuracies and their relevance 

to geologic mapping 

by 

L.U. Bender, M.H. Podwysocki, L.C. Rowan, and J.W. Salisbury 

U.S. Geological Survey 
Reston, Va. 22092 


1 

S 

l 


Landsat-4 Thematic Mapper (TM) scenes for Washington, D.C. 40109-15140, 
November 2, 1982) and Macon, Georgia (40050-15333, September 4, 1982) have 
been examined to determine their radiometry and geometric accuracy. In 
addition, these two scenes also are being analyzed to determine the ability to 
identify specific rock types with the added near Infrared TM bands located at 
1.6 and 2*2 urn. 

Each of the seven TM bands of the geometric* lly corrected (P-tape) 
Washington, D.C. scene were contrast stretched and recorded on film using an 
Optronics P-1500 film writer. Thirty-eight control points were visually 
identified from 1 :24, 000-scale topographic maps and a 40-micrometer circular 
hole was drilled at each control point on the Images of bands 1, 3, 5 and 7. 
Subsequently, rectangular coordinates of each control point were measured with 
a photogramme trie comparator with a resolution of 1 micrometer. Two 
transformations, the similarity and the affine, were fit by the least squares 
method to the comparator coordinates and to their respective UTM coordinates 
as taken from the maps. The similarity transformation only allows the data 
points to be translated, scaled and/or rotated about the coordinate system 
whereas the affine translation also allows for separate linear scaling of each 
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coordinate axis. Table 1 shows the vector residuals resulting from each 
transformation for several of the bands. 


Table 1 Root Mean Square Error of Vector Residuals 


Band 

Similarity 

Affine 

1 

72 m 

45 m 

3 

70 m 

45 m 

5 

68 m 

48 m 

7 

69 m 

45 m 


It Is apparent that the application of the affine transformation 
significantly reduces the root mean square error In all cases. One possible 
cause may be the film writing device. A highly significant skewness was 
observed In previous tests of grids printed from the same device as used In 
this study. This was possibly caused by lack of parallelism between the axis 
of the rotary drum and the axis of travel of the light source. 

Band to band registration was checked for the Washington, D.C. scene 
viewing several small portions of the scene on an Interactive video display 
device. Several strongly contrasting targets, such as airport runways and 
small ponds, were used to visually determine the offsets. TM bands 1-4 were 
perfectly coregistered. Using channels 1-4 as references, bands 5 and 7 were 
found to be displaced 1-2 pixels eastward. Band 6 was displaced 4 pixels 
eastward, which effectively means It Is displaced by one of Its 120m 
resolution units. 
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A geometrically uncorrected version (A-tape) of the Macon, Georgia scene 


was examined to determine some of the radiometric characteristics of the 
data. This Included analysis of the 16 Individual detectors In bands 1-5 and 
7 and each of two detectors of channel 6. Systematic gaps were found within 
the limited range of digital numbers recorded for the scene. These can be 
attributed to the integer truncation of decimal numbers Involved In the 
calibration procedure during production of the computer compatible tapes. 

The potential of sensor striping caused by the separate gains and offsets 
required for each of the 16 sensors of the visible and near lnfared TM bands 
was investigated using the A-tape of the Macon scene. The two sensor thermal 
Infrared band also was examined. Means and standard deviations were 
calculated on a modulo-16 basis. The results are summarized In table 2. The 
range of means is typically less than that of Landsat MSS values when the MSS 
data are scaled into 8-bit space. Contrast stretched images of portions of 
the Macon scene were recorded onto film and visually examined for evidence of 
striping. Striping is not evident In channels having sensors with a range of 
mean differences less than 0.4 between their high and low values. Striping 
becomee apparent when difference values reach 0.5 and becomes objectionable 
with values greater than 0.6. The means and standard deviations of TM band 1 
appear excessively high In comparison with TM channels 2 and 3, suggesting 
that gain and/or offset values for these sensors may be set too high. It Is 
unlikely that these Increased values in TM channel 1 can be accounted for by 
atmospheric scattering alone. 

A modulo 32 striping also was noted for some of the TM bands over 
homogeneous targets such as water bodies. The striping appears clustered in 
two groups of 16 adjacent scan lines and probably is related to the alternate 
back and forth sweeps of the pushbroom scanner. Additional analysis is 
undervay. 
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Table 2-Comparison of Digital Number Statlatlca for Modulo-16 Striping 
of TM data 


TM Channel 

1 

2 

_3 

4 

_5 

6 

_7 








range of 

85.5 

34.2 

30.6 

74.5 

66.8 

126.0 

22.7 

sensor means 

86.4 

34.8 

31.2 

74.9 

67.1 

126.7 

23.1 

dif ference 

0.9 

0.6 

0.5 

0.4 

0.3 

0.7 

0.4 

range of 

standa rd 

17.649 

9.067 

11.295 

16.665 

23.449 

15.527 

11.679 

deviations 

18.076 

9.504 

11.661 

17.026 

23.931 

15.980 

12.271 

dif ierence 

0.427 

0.437 

0.366 

0.361 

0.482 

0.453 

0.592 

median 

sensor value 
standard 

86.0 

34.6 

30.95 

74.6 

67.0 

126.7 

22.9 

deviation 

17.85 

9.24 

11.468 

16.897 

23.565 

15.963 

11.979 


Band ratlolng of specific TM bands was done for both the Macon, Georgia 
and Washington, D-C. scenes. Specific band ratios were used to determine 
their effectiveness in defining vegetation and rock mineralogy. The TM 3/4 
band ratio was selected to define vegetation and the 5/2 band ratio was used 
to define limonitic rocks. The 5/7 band ratio was evaluated to test the 
effectiveness of the new 1.6-and 2.2-^im bands in detecting clay and carbonate 
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bearing rocks. Preliminary results show that clay minerals containing 
absorption bands at 2.2 pm, which are centered In TM band 7, can be readily 
detected when the band ratio linages are suitably contrast stretched. Bright 
rocks containing carbonates, such as marbles, which contain absorption bands 
In the 2.33-/im region, near the long wavelength cutoff of TM band 7, are only 
evident upon harsher contrast stretching of the 5/7 band ratio. 

Geobotanical studies using the early November Washington, D.C. scene show 
that monospecific forest canopies can be differentiated. A false-color 
composite Image of TM bards 1,4 and 7 allowed the photointerpreter to 
distinguish between relatively pure sta.ids of Virginia pine (Plnus vlrginica ), 
loblolly pine (£. taeda ) and white pine (?_. strobus ) , although de-'se stands of 
Virginia cedar ( Juniperus vlrginica ) could not be distinguished from the 
Virginia pine. All groups of coniferous trees could be distinguished from 
deciduous canopies. 


5 





JSC THEMATIC MAPPER SCROUNGE DATA PREPROCESSOR 


J. Gilbert 

NASA Johnson Space Center 
Houston, Texas 


ABSTRACT 

With the launch of Lemdsat-D In July 1982, the remote sensing community 
has been provided a near Earth-orbiting multlspectral sensor which has 
Increased spatial resolution and superior spectral discrimination than 
any previously launched Landsit Instrument. This new Instrument, the 
Thematic Mapper, Is a forward step In the progression of remote sensing 
technology. Incorporating the knowledge and experience gleaned from early 
scanning Instrument endeavors as well as utilizing the latest available 
satellite stability and pointing capability. 

The Immediate concern of Landsat ground processing segments became the 
accommodation of a data volume output from the Thematic Mapper which Is 
seven times as great as multlspectral scanner output for the same scene 
coverage. Often, existing multlspectral scanner processors were In-place 
since the early 1970 ' s and were not capable of handling such volumes of 
data. The problem consequently became a manifestation of the adaptation of 
existing systems, wherever feasible, and In many Instances, the Implemen- 
tation of new processors for Increased data volume handling. 

This paper presents the development of the Thematic Mapper ground pre- 
processing capability of the Earth Observations Data Laboratory (EODL) at 
Johnson Space Center. The Implementation of new processing elements and 
the utilization of existing systems In the EODL to handle Thematic Mapper 
data are reviewed in detail. The EODL Thematic Mapper preprocessor has been 
designed to provide AgRISTARS researchers with satellite Imagery data over 
selected areas of interest In data sets which are sized to be computationally 
manageable. In addition, the EODL Thematic Mapper preprocessor Includes 
the capability to provide Image products In support of research activity. 



Processing of Thematic Mapper Simulator Data 


Stephen Ungar 

NASA/Goddard Space Flight Center 
Code 923 

With the advent of advanced satellite-borne scanner systems, the geometric 
and radiometric correction of aircraft scanner data has become increasingly 
important. These corrections are needed to reliably simulate observations 
obtained by such systems for purposes of evaluation. This paper reviews two 
approaches to reducing distortion in aircraft scanner dati: by a combination 

of local and global "rubber sheet" polynomial transformations based on ground 
control points; and by modelling the effect of aircraft motion on the scanner 


scene. 




ABSTRACT 

Evaluation of Land sat-4 Thematic 
Mapper Ground Segment Par f oraanoe 
Without Ground Control* 

N. A. Bryant, A. L. Zobrlst, P. C. Billingsley 
S. Z. Frledean, B. Gokhman, T. L. Logan 

An analysis was performed on the P-format oosputer oospatlble tapes for the 
Landsat-4 Thesatlc Mapper (TM) scene 40109-15140 acquired November 2, 1982 
over Washington, D. C. and environs. Three tests of sensor geometry were 
undertaken: a) band-to-band registration, b) line-to-line registration 

between swaths, and c) geometric correction of sensor and spacecraft/epheaerls 
characteristics. The band-to-band and line-to-line registration was measured 
at one hundred pixel s pacings along a line using the phase correlation image 
alignment method developed by Kuglln and Hines and adapted to a one- 
dimensional FPT correlation technique. TLe sensor and spacecraft geometric 
calibration analysis was checked by identifying 75 ground control points in 
the TM scene and on 1:24,000 topographic maps. A least squares fit of the 
gcps was computed and the vector offsets of the residuals was tabulated and 
plotted. 

An analysis of the tabulated and plotted results revealed the following 
characteristics: a) Band-to-band registration in the primary foeal plane 

(bands 1-4) is less than 0.1 pixel with the possible exception of band 4, 
which nay be slightly above 0.1 pixel, b) Band-to-band registration in the 
secondary focal plane (bands 5 A 7) is less than 0.1 pixel, c) Band-to-band 
registration between bands in the primary and seoondary focal planes varies 
between minus 0.1 and minus 0.35 pixels on the average where good correlations 
oould be obtained (note the different spectral responses for each band made a 
comparison difficult sometimes), d) Line-to-line mis-reglstratlon averages 
less than 0.1 pixel, with some apparent local mis- registration of 0.3 pixel. 
Plrts failed to show any systematic mis-registratlon effects that can be 
directly associated with local jitter, e) Ground oontrol point offsets from 
the residual surface fit are generally good (RMSE of 4.43 pixels or 132.77 
meterj), although some points are off to a significant degree. An analyses of 
the TM P-data shows that a slow spacecraft roll oo edition of under one hertz 
exists in this scene. It is expected that an even distribution of gcps oould 
effectively remove the roll distortion encountered. 

•This abstract presents the results of one phase of research carried out at 
the Jet Propulsion Laboratory, California Institute of Technology, under 
contract No. NAS7-100, sponsored by the National Aeronautics and Space 
Administration. 
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TM RESULTS SESSION 

EARLY EVALUATION OF TM GEOMETRIC CORRECTION SYSTEM 


Eric Beyer 
General Electric Co. 


Results of available Thematic Mapper Geometric Correction System evaluations 
are given. The early mission attitude control performance Is described. 

This Includes attitude deviations which result from Thematic Mapper, Multi- 
spectral Scanner, solar array and TDRSS antenna operations. Thematic Map- 
per performance is given with respect to active line length, total scan 
time and scan nonlinearity. As ground control points are not yet used in 
processing, evaluations of the processed Imagery Is limited to forward-to- 
reverse scan continuity and band-to-band registration. 



The Use of Linear Feature Detection to Investigate Thematic 
Mapper Data Performance and Processing 


Dr. C. Gurney, Mr. K. Eng 

Systems and Applied Sciences Corporation 

5809 Annapolis Road 

Hyattsville, Maryland 20784 


Abstract 

Linear features are defined as long, narrow features such as roads, rivers and bridges. In 
addition to applications utility for mapping, they may be used to investigate the 
geometric properties of remotely sensed data. 

Data suitability for feature detection is investigated by considering the contrast between 
typical feature and background materials. 

Applications of the detector model at a series of thresholds, and estimation of the 
resulting detection accuracies for features of known width and composition, allows 
estimation of the ground IFOV (pixel size). 

Band to band registration is estimated by considering the relative positions of features 
between bands. Results indicate an off-set of about .75 pixels along scan and .5 pixels 
across scan between features detected in the primary focal plane (bands 1-4) and those 
in the cold focal plane (bands 5 - 7). 
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Spatial Resolution Estimation of Landsat-9 
TM and MSS Data 


By 
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SPATIAL RESOLUTION ESTIMATION OF LANDSAT-4 
TM and MSS DATA 


Clare D. McGlllem 
Paul E. Anuta 
Kai-Bor Yu* 


Laboratory for Applications of Remote Sensing 

and 

School of Electrical Engineering 
Purdue University 
West Lafayette, IN 47906 


SlftMARY 

In order to verify that the Landsat-4 sensors are operating within 
specifications. It is useful to estimate the system parameters by analysis 
of the measured data. One parameter of particular Interest Is the sensor 
point-spread function (PSP) which determines the resolution of the system. 

A method of estimating the PSF has been developed that utilizes data 
obtained during scanning of ground elements having Identifiable geometric 
and radiometric structure. These data are then processed In such a manner 
as to recover either the PSF Itself or to estimate the parameters of an 
assumed functional representation of the PSF. 

The measured data can be expressed In the spatial domain as a convolu- 
tion of the scene with an overall point-spread function: 

g(x,y) - h(x,y) * f(x,y) 


* 


Currently at Virginia Polytechnic Institute, Department of Electrical 
Engineering, Blacksburg, VA. 
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h(x,y) Is the overall point-spread 

funcf ' -m of the sensor system 

g(x,y) 1 a the resulting image 

Given g(x,y), we wish to determine h(x,y) or its Fourier transform H(u,v). 
To do this, some deterministic element of the input f(x,y) must be known or 
assumed. Although the theory takes into account the two-dimensional nature 
of the element, the Initial experiments have been limited to the one-dimen- 
slonal case. Three useful elements of this type are: 

1. A step function represented by an abrupt change in gray level along a 
row or column of the data. 

2. An Impulse represented by a narrow-width discontinuity along a row or 
column of the data. 

3. A rectangular pulse represented by a sequence of two steps in opposite 
directions along a row or column of the data. 

The point-spread function estimation method described in the paper 

uses horizontal and vertical roads in Landsat-4 MSS and TM images. The two- 

dimensional function is solved as two one-dimensional functions under the 

assumption that the PSF is separable; i.e. , 


h(x,y) - h (x) h (y) 

* / 

The steps in carrying out the solution are: 

1. Each scan line across a road is modeled as an ideal point response 

f i P i " P i 6(x-x i ) 

P i - magnitude of the ideal point source 
x^ ■ location of the point source 


2. Each point response is adjusted such that the locations of the point 

sources align with each other and an average point response is computed. 
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8l (x) - p A 5(x-x 1 ) * h(x) - p A h(x-x 1 ) 
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M MM 

g(x) -is g (x+x ) • 15 E P 4 h(x) - kj E P ] h(x) 
" i -1 1 1 M 1-1 1 M i -1 1 


Therefore : 
h(x) - K g(x) 

1 M 

where K - l/[± E P.] 


Once an estimate for the PSF has been obtained to within an (unknown) 
constant factor K, a quantification of the width of the function can be 
made. Three methods are of general Interest and are defined In one dimen- 
sion, as follows: 

1. Half amplitude width W, - width at which the magnitude of h(x) falls to 
one-half of its value 2t the origin. 

00 

_l h(x) dx 

2. Equivalent width W • 

* h 

max 

3. rms width W - 2 r 

rms 

00 

2 * 2 h(x) dx 

r • 

00 

h(x) dx 

Similar quantltltes can be defined to estimate the bandwidth of the 
modulation transfer function (MTF) and are represented by B^, B # , and 

The method was applied to Landsat-4 MSS and TM Images from the Webster County, 
Iowa area and the Chicago, IL area using interstate highways and major roads 


-4- 


ORlGtlNlAL PAGE IS 

Of. POOR QUALITY 


as test Input*. The reaults for these tests sre obtained In the form of 
both spatial widths of the estimated PSF in pixels and the bandwidth of the 
modulation transfer function in cycles per pixel and by plots of the estimated 
PSF and MTF. 

One result using 50 scan lines across a north-south road in Webster 
County, Iowa (TM data) is shown in Figure 1. The PSF width results are 
listed below: 

Vertical road in Webster Co. , Iowa area: 

W, - 2.0 

S 

W - 3.58 
e 


(] 

c 

$ 

1 


- 5.23 


Evaluation of the method using an assumed step response will also be 
presented, along with results from the MSS and using data both from Iowa 
and the Chicago areas. 


1 

1 

0 
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AN ANALYSIS OF THE HIGH FREQUENCY 
VIBRATIONS IN EARLY THEMATIC MAPPER SCENES 


John Kogut 
Eliane Larduinat 

Research and Data Systems, Inc. 
10300 Greenbelt Road 
Lanham, Maryland 20706 


The potential effects of high frequency vibrations on the final Thematic 
Mapper (TM) image have been evaluated for 26 scenes between day 4 and day 50 
after LANDSAT-4 launch. The angular displacements of the TM detectors from 
their nominal pointing directions as measured by the TM Angular Displacement 
Sensor (ADS) and the spacecraft Dry Rotor Inertial Reference Unit (DRIRU) give 
data on the along scan and cross scan high frequency vibrations present in each 
scan of a scene. These measurements were used to find the maximum overlap and 
underlap between sucessive scans, and to analyze the spectrum of the high 
frequency vibrations acting on the detectors. 

The history of the scan overlap and underlap for this time period shows that & 
consistent underlap which was present immediately after launch disappeared after 
day 10. The maximum overlap/underlap within a scene stabilized after day 15 from 
launch to an average of -0.337 pixels / +0.358 pixels for the reflective bands, and - 
0.034 pixels / +0.141 pixels for band 6. 

The Fourier spectrum of the along scan and cross scan vibrations for each 
scene was also evaluated. In the scenes analyzed, the along scan high frequency 
vibrations were strongly concentrated at 7 Hz and had amplitudes of about 80 
microradians. Additional vibration peaks were present in the 56 Hz to 77 Hz range. 
Cross scan vibrations were also concentrated at 7 Hz and had amplitudes near 7 
microradians. For the cross scan data additional spectral peaks occurred at 21 Hz 
and in the 56 Hz to 67 Hz range. Differences in structure and peak amplitudes 
between the spectra of the forward and reverse scans were also observed. 


Assessment of Thematic Mapper Band to Band Registration by the 
Block Correlation Method: Preliminary Results 

Don H. Card and Robert C. Wrigley 
NASA Ames Research Center 
Moffett Field, California 94035 

Frederick C. Mertz and Jeff R. Hall 
Technicolor Government Services 
Ames Research Center 
Moffett Field, California 94035 


A Thematic Mapper Arkansas scene (August 22, 1982) is being 
processed on Ames Research Center's block correlation program to 
assess band to band registration. This is a Fast Fourier 
Transform based program implemented on a CDC 7600 computer that 
correlates spectral values between bands using a 64 X 64 pixel 
block in the primary band image versus a moving 32 X 32 pixel 
block in the secondary band image. Correlations are performed on 
gradient images in order to minimize radiometric differences 
between bands. Initial results on band 3 versus band 5 indicate 
an approximate 1/2 pixel offset in the vertical direction and 1/4 
pixel offset in the horizontal direction, with standard errors of 
.3 pixel in each direction. Further results with other band 
comparisons will be presented. 
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INVESTIGATION OF TM BAND-TO-BAND REGISTRATION USING THE 
JSC REGISTRATION PROCESSOR 

S. S. Yao & M. L. Amis, l.ockheed-EMSCO 


SUMMARY: 

Several TM scenes from the Goddard "SCROUNGE" processing system were 
available at JSC for a quick-look Investigation. One topic of Interest is 
the accuracy of the TM band-to-band registration. It Is known that the 
detectors for the seven-band TM system are placed In two groups aboard 
the spacecraft. Located at the primary focal plane r.re detector arrays 
for band one to band four, while the remaining detectors are placed at 
the cold focal plane. Software In the SCROUNGE system Is used to register 
the bands together t*> an accuracy of less than one-fifth of a pixel, or 
six meters on the ground. It Is expected, therefore, that the data from 
the first four bands are well registered to each other. In the current 
Investigation, the JSC registration processor Is used to study the registra- 
tion between all TM bands, with special attention given to band 6, the 
thermal band, which has a four-times coarser resolution compared to the 
other bands. 

The JSC registration processor Is an Improved offspring of Goddard's 
Master Data Processor and the LACIE proceisor. It performs scene-to-scene 
(or band-to-band) correlation based on e..ge Images, which are derived 
from a percentage of the edge pixels calculated from the raw scene data, 
excluding clouds and other bad data In the scene. Correlations are perfor- 
med on patches of the edge Images, and the correlation peak offset In each 
patch Is estimated Iteratively to fractional pixel location accuracy. 

Peak offset locations from all patches over the scene are then considered 
together, and a variety of tests are made to weed out outllners and other 
Inconsistencies before a distortion model Is formed. Thus, the correlation 
peak offset locations In each patch will Indicate quantitatively how well 
the two TM bands registered to each other over that patch of scene data. 

When multi temporal acquisitions of TM data over the same ground track 
were made available, the JSC registration processor was also used to 
multi temporally register one acquisition to the other. Some multi temporal 
TM registration results together with the preliminary evaluation of single 
acquisition band-to-band registration results will be presented. 



GEODETIC ACCURACY OF LAND SAT- 4 MULTISPECTRA L SCANNER AND THEMATIC MAPPER DATA 
Principal Investigator: June M. Thornodagard 

EROS Data Center , U.S. Geological Survey 

The EROS Data Center, as the primary facility for generation and dissemination 
of Landsat data products, must evaluate the geodetic accuracy of data from both 
the Multlspectral Scanner (MSS) and Thematic Mapper (TM) processing systems. 

The standard Landsat-4 MSS and TM digital tapes contain information which 
relates the map projection coordinate system to the image coordinate system. 
Therefore, by converting the latitude and longitude of a given location to the 
map projection coordinates, the image coordinate location can be calculated. 
Verification of the Image's geodetic accuracy can then be accomplished by 
comparing the calculated image coordinates with the visually determined image 
location. 

Landsat-4 digital images produced by NASA are either system-corrected or 
precision-corrected (ground control points applied). This method of geodetic 
verification has been applied to images of both precision-corrected and 
system-corrected MSS and TM data. 






Thematic Mapper Quick-Look Analysis of the Washington, DC Scene (11/2/82) 

by 


James R. Irons 
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Darrel L. Williams 

NASA/Goddard Space Flight Center 
Earth Resources Branch/Code 923 
Greerbelt, Maryland 20771 


Richard S. Latty 
University of Maryland 
Department of Civil Engineering 
College Park, Maryland 


Mark L. Stauffer 

Computer Sciences Corporation 

Silver Srings, Maryland 


ABSTRACT 

Preliminary analysis of the landsat-4 Thematic Mapper (TM) scene of Washington, DC 
(November 2, 1982) has been completed. The objective of this effort » as to 
provide very near term quantitative results on assessment of the increased 
information content provided by the improved spatial, spectral, and radiometric 
resolution of the TM relative to the Multispectral Scanner (MSS). 

Four different “data sets" were analyzed. These consisted of: (a) the original 

TM data (28.5 meter pixel, 8 bits, 6 channels), (b) the original data degraded 
to 6 oit quantization, (c) a 3 channel subset of the original data (TM bands 
2, 3, and 4 which are similar to MSS bands 1, 2 and 4), and (d) the original 
data degraded to MSS spatial resolution. [The thermal data (band 6) were 
excluded from this study due to the difference in spatial resolution (120 
meters.)] Nine different study sites, 25b by 256 pixels in size, were randomly 
selected from the full scene for analysis. Recent aerial photography was 
available and all sites were visited during the last week of October to collect 
ground reference data. All nine study areas for each of the four data sets 
were classified using a per-point guassian maximum likelihood classifier. 

The following results were obtained: (a) the reduction of quantization level 

from 8 bit to 6 bit caused a decrease in overall accuracy ( 7%), (b) the use 
of only three bands (TM 2, 3 and 4) covering the visible and near infrared 
portion of the spectrum, caused a decrease in overall accuracy ( 7%), and (c) 
the decrease in spatial resolution res.lted in an increase in overall accuracy 
( 4X). Results (a) and (b) indicate that the increased radiometric and spectral 
resolution of the TM instrument do provide increased information content. 

The results of the spatial resolution degradation is somewhat misleading, in 
that the result is more a function of the type of classifier used (i.e., per- 
point), rather than a function of spatial resolution. This result points to 
the need for new classifiers, such as contextual classifiers, which take into 
account the increased spectral heterogeneity in higher resolution data. 
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Preliminary Comparisons of the Information Content and Utility 

of TM Versus MSS Datu 


by 

Brian L. Markham 
Earth Resources Branch 
NASA/GSFC 

Greenbelt, MD 20771 
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Comparisons were made between comparable subscenes from the first TM scene 
acquired of the Washington, DC area (July 29, 1982-4 bands only) and a MSS 
scene acquired approximately one year earlier (July 11, 1981). The Landsat-4 
MSS scene of July 29, 1982 was never made available. Three types of analyses 
were conducted: a water body analysis, a principal components analysis and a 

spectral clustering analysis. 

The water body analysis compared the capability of the TM to the MSS for detecting 
small uniform targets. All water bodies (>10 meters) were located on aerial 
photographs collected July 13, 1982 and categorized by the size of the maximum 
inscribable circle. Each scene was independently clustered and the clusters 
assigned to water or non-water. Each water body was said to bo detected for 
the MSS or the TM if a minimum of one pixel at its location was classified 
into a water cluster. Of the 59 ponds located on the aerial photographs 34 
(58%) were detected by the TM with six commission errors (15%) and 13 (22%) 
were detected by the MSS with three commission errors (19%). The smallest 
water body 'detected* by the TM was 16 meters; the smallest 'detected' by the 
MSS was 40 meters. 

For the principal components analysis, means and covariance matrices were 
calculated for each subscene, and principal components images generated and 
characterized. The 4-band TM image had slightly higher dimensionality than 
the MSS in this rural -suburban area, with 2.6% of the variance in the third 
component as compared to 1.2%. The first two principal components of the TM 
and MSS were similar, containing brightness-like and greeness-1 ike features. 

The MSS third and fourth principal components were dominated by noise and 
striping whereas the third TM principal component (a red to blue-green contrast) 
was less noisy and appeared to contain information useful for separating built- 
up features from bare soil. This third component detected the spectral flatness 
of construction materials in the visible spectrum as opposed to the sloping 
response of bare soils. 

In the spectral clustering comparison each scene was independently clustered 
and the clusters were assigned to informational classes. The clusters were 
assignable to the same classes on the 4-band TM as on the MSS data with the 
exception that the TM data provided sepa» 'tion of bare soil areas from areas 
of buildings that MSS data did not provide. The TM band 1 contributed to this 
improved separation. 

The preliminary comparisons indicated that TM data provides enhancements over 
MSS in terms of (1) small target detection and (2) data dimensionality (even 
with 4-band data). The extra dimensior partially resultant from TM band 1, 
appears useful for bui lt-up/non-bui It area separation. 
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INFORMATION CONTENT OF DATA FROM THE LANDSAT 4 THEMATIC 
MAPPER AND MULTISPECTRAL SCANNER 
John C. Price 

USDA, Hydrology Laboratory 
Beltsville Agricultural Research Center 
Beltsville, Maryland 20705 

The word "information" has both a popular and a technical meaning, the 
latter referring to the statistical variability associated with a data set. 
This study addresses the question of the quantitative increase in the 
information content (in the technical sense) of Land "tat 4 Thematic Mapper (TM) 
data as compared to data from the Multispectral Scanner (MSS). 

An intuitive concept of the information content of image data may be 
gained from several examples: 1) A uniform featureless scene, such as a desert 
area, has a low information content while a "busy" scene, such as a city, 
contains more information; 2) A scene having no spectral variability, i.e. all 
wavelength channels having a common brightness value at each point - a black 
and white scene, has a lower information content than one in which brightness 
values differ from one spectral channel to another in a spatially varying way 
- a multicolored scene; 3) A data set obtained at low spatial resolution has 
many features blurred or averaged out, and thus has a lower information 
content than the equivalent scene observed at higher resolution. 

In the present case, the TM data contain more information than the MSS 
data because of 1) the greater number of spectral channels, 7 versus 4, 2) the 
improved spatial resolution, 30m versus 82m, and 3) the increased precision 
resulting from better detector response and decreased quantization intervals, 
e.g., 8 bit versus 6 bit data. 
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The analysis procedure consists of the application of standard mathematical 
methods (principal components, autocorrelation, delta transformat ion, etc.) to 
data sets from the MSS and TM. The subject area, information theory, as well 
developed and software already has been developed to support the study. 

Computer runs have been carried out for several subareas from a data set 
acquired simultaneously by TM and MSS over a test area in northeast Arkansas. 
(Scene center North 36°Q2 ' , West 90°Q2 ' , date, August 22, 1982) This 
scene is one of th- most agriculturally diverse areas in the U.S. As one 
would expect, the 6 visible-near IR channels of the thematic mapper provide 
more information than the 4 channels of the mul tispectral scanner. A rough 
estimate of 20 bits per pixel for TM, 10 bits per pixel for MSS has been 
computed for these subareas. The*e numbers will be revised downward when 
allowance is made for noise in the data. 

Since the 6 reflective channels of the TM have a capacity of 48 bits per 
pixel (6 channels x 8 bits per channel) while the 4 MSS channels have a 
capacity of 27 bits per pixel (3x7 bits ♦ 1 x 6 bits) one can readily 
compare the information gathering efficiency of the two instruments. It 
appears that the ratios (information per pixel)/(possible information per 
pixel) are approximately equivalent - 2 0/48^10/2 7, with noise not yet 
accounted for. If data utility is related to the quantitative estimate of 
information then the greater spatial and spectral capability of the TM 
represents a substantial improvement in remote sensing technology. However, 
the present assessment is preliminary and study of additional data sets is 
required. 
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LANDSAT 4 INVESTIGATION OF THEMATIC MAPPER 
AND MULTISPECTRAL SCANNER APPLICATIONS 


By D. T. Lauer- and F. A. Walta— 


ABSTRACT 

Tha objective of this experiment was to determine the Incremental value 
of Landsat 4 thematic mapper (TM) data over multlspectral scanner (MSS) data 
for acquiring resource Information from selected test sites. Only data 
collected on November 2, 1982, over the Washington, D.C., test site have been 
analysed. 

Data from the six visible bands (1, 2, 3) and near-infrared spectral 
bands (4, 3, 7) of the TM data were used for several land use applications. 

The statistical characteristics of the Information obtained from the various 
bands were determined by principal component analysis, correlation analysis, 
and feature classification. Preliminary results of these analyses showed that 
the most Informative combination of bands Includes the near-infrared band 4 
with at least least two other bands, either visible or Infrared. This 
Information was used to select bands to be uaed In a hue-intensity, 
saturation-image transformation algorithm. 

A second approach based on the relative lnterpretablllty of features In 
the black-and-white products was used to make optimum band selections for 
falae-color Image products. 

Work was begun building a multilayered digital data base covering the 
Washington, D.C. , area. TM and MSS data were registered with digital 
elevation data (digital terrain tape), and TM conjugate stereoscopic Images 
were then made by Introducing parallax via the digital elevation data. Into 
all seven bands. Both true and false color stereoscopic Image products were 
made for the Washington, D.C., test site. Also, the digital data were uaed to 
make three-dimensional perspective views In which the axlmuth (viewing angle) 
and location (height) of the obeerver were varied. 

Digital products were statistically analysed and photographic products 
manually Interpreted by experienced EROS Data Center scientists to determine 
the Incremental value of Landsat 4 TM over MSS data. This value was then 
related to the physical parameters of the two systems. 


— U.S. Geological Survey 

2 / 

— Technicolor Government Services, Inc. Work performed under 
U.S. Geological Survey contract 14-08-0001-20129. 
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PRELIMINARY EVALUATION OF THEMATIC MAPPER 
DATA FOR RENEWABLE RESOURCES 


D. E. Pitts 

NASA Johnson Space Center 
Houston, Texas 

J. Carnes, E. Reyna, R. B. Pollock, and C. Sorensen 
Lockheed Engineering and Management Services Company, Inc. 

Houston, Texas 

ABSTRACT 

The Thematic Mapper data are being used to empirically determine the 
effect of atmospheric scattering on spatial resolution. The response 
of the Thematic Mapper to discontinuous changes In scene reflectance 
was used to estimate the system optical transfer function. These 
estimates are being compared with protofllght measurements to determine 
If the observed spatial response Is consistent with atmospheric 
scattering model predictions of spatial resolution. 

In order to evaluate the effect of increased resolution on classification, 
an unsupervised clustering algorithm (CLASSY) was applied to the TM data 
and to simultaneously acquired MSS data. By observing features which were 
not resolved Into separate clusters using TM data, an Indication of how 
the Increased TM resolution would affect classification was obtained. 

The results of this evaluation Indicated that using TM data would produce 
separate classes for small features such as roads and homesteads, while 
using MSS data would not separate these sma 1 ! features. 
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THEMATIC MAPPER DATA QUALITY AND PERFORMANCE 
ASSESSMENT IN RENEWABLE RESOURCES/ AGRICULTURE 
REMOTE SENSING 

Robert M. Blzzell and Harold L. Prior 

NASA-Lyndon B. Johnson Space Center 
Houston, Texas 77058 


For the past decade the Johnson Space Center has been involved In the research 
and development of technology that utilizes satellite based remotely sensed data 
to perform global crop inventories. The Large Area Crop Inventory Experiment 
(LACIE) established the proof-of-concept with the LANDSAT MSS. The fundamental 
understanding gained with the LACIE formed the basis for the development of 
these techniques in the framework of a systems approach that has been evaluated 
in the Agriculture and Resources Inventory Surveys Through Aerospace Remote 
Sensing (AgRISTARS) project. The successes with these Joint projeots have been 
manifest by the Installation and periodic upgrade by the Foreign Agriculture 
Service (FAS) of the USDA within their ongoing operational forecasting system. 

The limits reached with the MSS emphasized the need for the Thematic Mapper (TM) 
to further expand the progress made to date. It is believed that the Increased 
spatial resolution will provide solutions to proportion estimation error due to 
mixed pixels, and the increased spectral resolution will provide for the 
identification of Important agricultural features such as crop stage, and 
condition. 

This paper will describe the results of analyses conducted relative to these 
hypothesis from sample segments (approximately 30 sq. miles) extracted from the 
4-band Detroit scene and the 7-band Mississippi, Ar engineering test scene. 
Several studies were conducted to evaluate the geometric and radiometric 
performance of the TM to determine data viability for the more pertinent 
investigations of TM utility. In most cases this requirement was more than 
sufficiently satisfied. This allowed the opportunity to take advantage of 
detailed ground observations for several of the sample segments to assess class 
separability and detection of other important features with TM. 

The results presented regarding these TM characteristics show that not only is 
the increased definition of the within scene variance captured by the Increased 
spatial and spectral resolution, but that the mid- IR bands (5 and 7) are 
necessary for optimum crop type classification. Both qualitative and quanti- 
tative results are presented that describe the improvements gained with the TM 
both relative to the MSS and on its own merit. 
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REMOTE SENSING INVESTIGATIONS OF WETLAND BIOMASS AND 
PRODUCTIVITY FOR GLOBAL BIOSYSTEMS RESEARCH 
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Michael Hardlsky 
Vytautas K1 etnas 
College of Marine Studies 
University of Delaware 
Newark, DE 19711 
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ABSTRACT 

Wetland ecosystems constitute a minor portion of the earth In an ^ 

areal sense, yet they are probably one of the most Important ecosystems 
when considering the global cycling of some elements having major reserves 
In the atmosphere. Among the outstanding characteristics of wetlands 
are the anaerobic sediments associated with these ereas. Within the 
anaerobic zone, fermentation, denitrification, sulfate reduction and 
methanogenesls provide the mechanisms for C, N and 5 to be returned to 
the atmospheric reservoir (Odum 1979). Considering that wetlands 
ecosystems are among the most productive (per unit area) on earth and 
that they maintain an environment favorable to the nutrient transformations 
necessary to return C, N and S compounds to the atmospheric reservoir. 

It Is critical that wetland ecosystems be evaluated as extremely Important 
areas In global ecology models. 


Systems analysis on a global scale will necessitate Information 
gathering from orbiting platforms. The vegetation of wetlands will be 
the most measurable surface feature. Since plants are good Indicators 
of edaphlc conditions, they can serve as monitors of the anaerobic 
conditions In the soil. Therefore, monitoring biomass of wetlands 
ecosystems can provide Information on net primary production and on the • 

chemical and physical status of wetland soils relative to anaerobic 
microbial transformation of key elements. 

Multlspectral remote sensing techniques have been successfully used 
to estimate r.ucrophytlc biomass in wetlands systems (Jensen 1980, Hardlsky 
et al. 1983a). The Improved spatial, radiometric and spectral resolution 
of the Landsat-4 Thematic Mapper over the Landsat MSS will greatly 
enhance multlspectral techniques for estimating wetlands biomass over 
large areas. These techniques can provide the biomass data necessary 
for global ecology studies. 


Approach 

A data base of spectral and biomass information has been compiled for 
the dominant salt marsh plant Spartlna alterni flora . Spectral data were 
collected with a hand-held GSFC Mark II radiometer (Tucker et al . 1981). 
The radiometer was configured with three spectral bands deemed most 
useful for studying plant canopies. A red band (0.63-0.69 yin) sensitive 
to chlorophyll density, a near Infrared band (0.76-0.90 ym) sensitive to 
tissue structure or biomass and a middle Infrared band (1.55-1.75 ym) 
sensitive to leaf moisture were used. These spectral bands correspond 
to bands 3, 4 and 5, respectively, of the Thematic Mapper. 
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This Initial data acquisition was designed to accomplish the 
following objectives: 

1. To determine the range of biomass density for naturally 
occurring S. alternl flora which can be discriminated using 
spectral data. 

2. To evaluate normal seasonal changes In several S. alternl flora 
stands with respect to spectral discrimination of these changes. 

3. To develop regression models equating spectral radiance and 
aboveground biomass of alternl flora for use as predictive 
models In future nondestructive estimates of biomass. 


If nondestructive biomass estimation of tidal wetlands using 
remote sensing Is to be a viable operational methodology, many of the 
plant species and plant associations other than S^_ al terniflora must 
be Investigated. By developing regression models for estimating biomass 
of the major salt marsh plant species and for the major plant associations 
o* typical brackish marsh systems, realistic estimates of total marsh 
biomass and production can be derived from remotely sensed spectral data. 

The central theme of future research Is to develop the methodology 
necessary for spectral estimation of marsh plant biomass and to use this 
methodology for the estimation of net primary productivity of tidal 
wetlands. Both ground-based and airborne platforms for spectral data 
collection will be tested and marshes experiencing different salinity 
regimes will be Investigated. The principal goals of this research can 
be summarized as follows: 

1. To develop predictive regression models equating spectral 
radiance and plant biomass for a variety of marsh plants 
exhibiting different leaf morphology, growth form and phenology. 

2. To employ the regression models to estimate aboveground plant 
biomass temporally for salt and brackish marshes and to estimate, 
solely from radiance data, the net primary production of the 
marshes. 

3. To compare biomass estimated using spectral data collected 
from the ground, from low altitude aircraft and from Landsat 
MSS and Thematic Mapper with biomass estimated using harvest 
techniques. 

4. To expand the spectral estimation technique to coastal mangrove 
swamps . 
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Results 


Results from earlier studies suggest that using spectral data, one 
can successfully detect changes In S. alternlflora biomass (Hardlsky et 
al. 1983b, c). In addition to monitoring perturbation, regression 
models have also been developed for predicting S_ alternlflora biomass 
over an entire growing season from spectral dita (Hardlsky et al . 1983d). 
Models developed over the growing season Include seasonal variations In 
biomass density and Illumination Intensity. 

Using the models, we then collected an Independent set of biomass 
and spectral data for alternlflora . We estimated standing crop 
biomass and net primary productivity using the regression model and then 
compared the remote sensing estimates to the harvest estimates of biomass 
(Hardlsky et al. 1983a). The results suggest biomass and net primary 
productivity estimates for S. alternlflora marshes using multlspectral 
data are possible. Albeit these initial Investigations have relied 
solely on spectral data collected from the ground, the spectral regions 
and algorithms suggest that successful estimation of wetlands biomass 
and production from the Thematic Mapper will be possible. 

Preliminary data are now being collected for wetland plant species 
other than al ternlflora . Initial trends suggest that plants with 
different morphologies ani canopy structure yield a different association 
between spectral data and biomass. For example, gramineous canopies 
(like S. alternlflora ), leafless stem canopies (like Salicornia sp.) and 
broadleaf canopies (like Iva frutescens or Sol 1 da go sp.) exhibit differing 
reflectance per unit leaf biomass density. Canopy structure (orientation 
of leaves) plays an important role In determining spectral response. 
Adjustments to models must be made such that biomass estimations are 
possible for the variety of plant types inhabiting wetland systems. 


Implications 

Based on spectral data similar to the Thematic Mapper and collected 
from the ground, we can suggest the following: 

1. The spectral configuration of the Thematic Mapper Is such that 
estimates of wetlands biomass and/or productivity can be made 
using spectral data. 

2. Wetland plant communities having different canopy configurations 
will require different biomass-spectral data algorithms. Thus 
plant species having similar canopy structure can probably be 
incorporated into a single model for biomass estimation. 

3. With algorithms for the major plant types, large area estimates 
of biomass and productivity will be possible from the Landsat- 
4 Thematic Mapper. 
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ASSESSMENT AND COMPARISON OF LANDSAT TM AND MSS 
DATA FOR DETECTING SUBMERGED VEGETATION 

Steven G. Ackleson 
Vytautas Klemas 

College of Marine Studies 
University of Delaware 
Newark, Delaware 19711 


ABSTRACT 


The objective of this research is to assess the spaclal , spectral, 
and radiometric resolution of LANDSAT TM and MSS data for detecting 
submerged grass beds in Chesapeake Bay. Preliminary results show the 
relative effectiveness of TM bands 1, 2, and 3 and MSS bands 4 and 5 for 
discriminating between submerged vegetation and mud or sand. Only the 
spectral characteristics of each band are considered. A single-scattering 
radiative transfer model was used to generate a spectral quality index 
(SQI) for shot-noise-l imited detectors. As SQI Increases, so does the 
appamnt contrast between target reflectance and background reflectance. 

TM band 3 (630nm - 690 nm) and MSS band 5 (600nm - 700nm) were 
found to best discriminate between vegetation and sand In turbid water. 
Corresponding SQI values representing a water depth of 0.5 meters were 
174.3 and 154.8 respectively. TM bands 1 and 2 and MSS band 4 were all 
of significantly lower quality, with SQI values of 69.4, 70.6 and 71.7, 
respectively. Similar results were obtained using mud as the background. 
The major difference was that all SQI /alues decreased. 

For each of the two situations, sand background and mud background, 
the optical properties of the water column appeared to have a minimal 
effect upon the relative usefulness of each band with respect to increasing 
water depth. The optimum band for detecting vegetation in shallow water 
remained optimum for deeper water. An example is presented In which the 
optimum band does change with increasing depth. 
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A Preliminary Evaluation of Thematic Mapper Data to 
Measure Suspended Sediments in Lake Chicot, Arkansas 

Frank F. Schisbe, Jerry C. Ritchie and Glen o. Boatwright 
USDA-ARS, Durant , Ok., Beltsvi lle,Kd. , Houston, Tx. 

A comparison was made between ground data collected 
from Lake Chicot, Arkansas and thematic mapper data 
collected on September 23» 1932. A preliminary analyses 
of limited data indicates that thematic mapper data may 
be useful in monitoring suspended sediments and chlorophyll 
in lake with high suspended sediment loads. Total suspended 
loads ranged from 166 to 508 ppm. TK band 3 appears 
to be most useful to measure suspended loads. Bands 
1,2, and 4 also have useful information on suspended 
sediments. Digital data from band 2 is correlated with 
chlorophyll concentration. The ratio of band 1 -band 2/ 
band 1 + band 2 is also correlated with chlorophyll 
concentration. Looking at water data only, bands 1,2, and 
3 appear to have similar information. Band 3 is also 
related to band 4* Bands 5*6, and 7 appear to have 
independent information content but it v:as not related 
to ary of the ground data we collected. 
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Sr: cur Reflectance from Thematic Mapper 

JEFF DOZIER 

University of California 
Santa Barbara 

Lend set -4 Early Results Sympoamm 
NASA Coddard Spaceflight Center 
February 22-24, 1063 

Calculations of snow reflectance in all 6 TM reflective bands (i.e. 1, 2, 3, 4, 5. and 
7), using Wiscombe and Warren's (I960) delta -Eddington model, show that snow 
reflectance in bands 4, 5. and 7 is sensitive to grain size An objective in cur investiga- 
tion is to interpret surface optical grain size of snov * r spectral extension of albedo. 
Our results so far are encouraging It also appears that the TM filters resemble a 
"square-wave" closely enough that we can assume a square-wave in our calculations. 
We calculated integrated band reflectance over the response functions, using sensor 
data supplied by Santa Barbara Research Center. Differences between integrating over 
the actual response functions and the equivalent square wave were negligible. 

Table 1 gives characteristics of the Thematic Mapper, and, for background infor- 
mation, the MultispccLral Scanner and NOAA Advanced Very High Resolution Radiome- 
ter. In the radiance columns of the table, the quantization errors and saturation radi- 
ances of the sensor bands are i . tipared with the solar constant, integrated through 
the sensor response functions Solar constant spectral distributions are from Thekae- 
kara (1970), adjusted to fit the integrated values of Hickey et al (1960) The last 
column in the table expresses the sensor saturation radiance as a percentage of the 
solar constant, integrated through the band response function While snow will still 
saturate in band 1. bands. 2 and 3 will saturate only for nearer-normal illumination 
angles and bands 4. 5, and 7 will not saturate. The saturation problem is not nearly as • 
severe as with the MSS, so the TM can be used to measure snow albedo and thus allow 
be sin-wide energy budget snowmelt calculations 

Table 2 shows calculations of integrated reflectance for snow over all reflective TM 
bands, and water and ice clouds with thickness of 1mm water equivalent over TM bands 
5 and 7. In the blue and green bands (1-2) snow reflectance is not sensitive to gram 
size, so measurements in these wavelengths will reveal the extent to wkdch snow albedo 
is degraded by contamination from atmospherical aerosols, dust, pine pollen, etc. In 
the red and near-infrared, snow reflectance is sensitive to grain size but not to contam- 
inants, so grain size estimates in these wavelengths can b used to spectrally extend 
albedo measurements. In bands 5 and 7, snow is dark ana clouds are bright. Initial 
frames with both snow and clouds from the western United States show that the TM can 
distinguish between them. 

References 

Hickey, J.R., L L Stowe, H. Jacobowitz, P. Pellegrino, R H. Maschoff, F. House, anc T.H. 
VonderHaar, 1960, Initial solar irradiance determinations from Nimbus 7 cavity 
radiometer measurements, Science, 208. 261-283. 

Thekaekara, M.P., ed.. 1970, The solar constant and the solar spectrum measured from 
a research aircraft, NASA TR-R-351 

Wiscombe, W.J., and S.G. Warren, 1980, A model for the spectral albedo of s:;-. r, 1, Pore 
snow, Journal of the Atrnosptirric Sciericts, 37, 27! 2-2733. 


PRELIMINARY EVALUATION OF TW FOR SOILS 


David R. Thompson and Keith E. Henderson 
NASA Johnson Space Center 
Houston, Texas 


Thematic Mapper data acquired over Mississippi Country, Arkansas, were 
examined for utility In separlng soil associations within general level 
alluvium deposited by the Mississippi River. The 1.58-1.79ym (Band 5) was 
found to separate the different soil associations farlly well when compared 
to the USDS-SCS general soil map. The thermal channel also appeared to 
provide Information at this level. A detailed soil survey was available 
at the field level along with ground observations of crop type, plant 
hlght, percent cover and growth stage. Soils within the fields ranged 
from unlfcrm to soils that occur at patches of sand that stand out 
strongly against the Intermingled areas of dark soil. Examination of 
the digital values of Individual TM bands at the field level indicates 
that the influence of the soil Is greater In TM than It was In MSS bands. 

TM appears to provide greater detail of whlchin field variability caused 
by soils than MSS and thus should provide Improved Information relating 
to crop and soil properties. However, this soil Influence may cause crop 
Identification classification procedures to have to account for the soil 
In their algorithms. 


FACIAL ANALYSIS OF THEMATIC MAPPER PRODUCTS 
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GEOLOGIC UTILITY OF LANDSAT-4 TH DATA: 
DEATH VALLEY, CALIFORNIA 

Michael Abrams, Anne Kahle, Alan Gillespie, 
James Conel, Harold Lang 
Jet Propulsion Laboratory 
Pasadena, CA 


The present LIDQA study Is being conducted In several phases to quantify 
the performance of the TM vis-a-vis various geological applications. These 
phases Include: (1) analyses of the geological utility of the data with 

resoect to the Increased spatial resolution and number of bands (compared 
tf whe MSS); (2) analysis of geometric accuracy; (3) analysis of radiometric 
performance of the TM scanner. 

This report presents some preliminary results of geological analysis of 
TM data for Death Valley, California, and observations on TM data quality. 
Geometrically and radlometrlcal ly corrected data (CCT-PT) have been 
examined for scene E-40124-17495, acquired November 17, 1982. Overall, the 
TM data appeared to contain a marked Increase In geologically useful 
Information; however, a number of instrumental or processing artifacts 
may well limit the ability of the geologist to fully extract this Information. 

A 27 by 18 km area was extracted from the scene for detailed analysis. 
Color-enhanced band composites, principal components composites, and color 
radio composites were produced from the TM data. Lithologic Interpretation 
maps were prepared and compared to detailed published geologic maps. More 
detail and delineation of alluvial fan units w < evident on the TM data; 
most of the sal t/sllt/evaporl te units In the Deuth Valley playa were 
apparent on the Images; many of the outcropping volcanic and sedimentary 
rock types were separable based on differences In their spectral response. 
However, substantial shadowing due to the combination of low sun angle and 
severe topographic relief limited the ability to Intepret information In 
the mountainous areas. High sun angle, summer scenes would greatly reduce 
this problem. 

The same region was extracted from Landsat-2 MSS data and also from aircraft 
multi spectral thermal data. Displayed as color composites and principal 
component Images, the TM data were superior to the MSS data In all respects. 
The multi spectral thermal data, however, allowed separation of alluvial 
units and lithologic units not seen on the TM data. This Indicates the 
complementary nature of the two data sets for providing lithologic 
separations based on sensing different physical properties. 

Two significant radiometric problems were observed In the TM data. The 
mean signal levels differ appreciably between the forward and reverse 
scan directions, possibly due to Inaccuracies In the radiometric calibrations. 
This effect manifests itself as banding In harshly contrast-enhanced images. 
The second problem Is most pronounced In ratio images and appears as 
striping, most severe at the edges of the scenes. The striping consists 
of paired lines of high-variance noise, with a 17-11 ne periodicity. They 



occur at the seams between the forward and reverse scan swaths, where 
extension scan lines have been Interpolated to fill data gaps. It appears 
that the Interpolation procedure used to process this scene was faulty. 

The TM data were registered to a topographic base map. Results of the 
registration procedure Indicate that the error In geometric accuracy of 
the CCT-PT data Is less than 2 pixels. A more rigorous test of accuracy 
Is planned. 


Preliminary Study of Information Extraction of Landsat TM Oata 
for a Suburban/Regional Test Site 
by 

David L. Toll 

NASA/Goddard Space Flight Center 

Earth Resources Branch/Code 923 ORIGINAL PAGE IS 

Greenbelt, Maryland 20771 £ p00R 0UAUTY 

Objectives 

The primary objective is to quantify the information extractable from TM data 
for a suburban/regional test site. Specific objectives include (1) quantifying 
information in processed TM data of land cover classes, and (2) assessing 
important TM sensor parameters affecting cover class discrimination. 

Approach 

Landsat TM (July 29, 1982 4 band; and November 2, 1982 7 band) and MSS 
(Landsat-2, July 11, 1981; and Landsast-4, November 2, 1982) P-data were 
registered scene to scene for a 15 km^ study area encompassing Beltsville- 
Laurel, MD. The MSS data was analyzed for comparisons to TM. In order to 
study the effects of improved TM sensor parameters--spatial resolution, quantization 
level, and spectral regions--on land use/land cover discrimination, the TM 
data was preprocessed. For both of the July and November TM scenes the data 
was degraded from 28.5m to 85m spatial resolution, compressed from 8 bits to 6 
bits quantization, and spectrally subset from 6 bands (4 bands in July) to 3 
bands (TM2, 0.52-0.60 *.m; TM3 0.63-0.69* m; and TM4 0. 76-0.90 /cm) . Thus for 
each of the two dates there are three factors (spatial resolution, quantization, 
and spectral regions) yielding a 2x2x2 combinatorial design of eight possible 
data arrangements (i.e., 6 bands, 28.5m, and 8 bits; 6 band, 28.5m, and 6 
bits; . . 3 bands, 85m, and 6 bits). To assist evaluations, CIR (1:40,000) 

of July 13, 1982, was photointerpreted to seven land use/cover classes: water, 

forest, agriculture, extractive, transportation, commercial /industrial , and 
residential. For cover class evaluations, pixels of TM and MSS were selected 
through a stratif ied-systematic sample in order to obtain independent and 
unbiased observations. Results fror.i correlations, histograms, two-dimensional 
spectral class plots, clustered maps, principal components, transformed divergence, 
and per-pixel classifications were evaluated. 

Findings 

Results from the TM sensor parameter study indicate a significant increase in 
classification accuracy and transformed divergence for the improved 8 bit 
quantization level and added spectral regions TM1 (0.45-0. 52* m) , TM5 (1.55- 
1.75* k m), and TM7 (2.08-2.35* m) . However, results from the spatial degradation 
procedure exhibited a significant improvement in classification accuracy and 
transformed divergence when using 85m vs 28.5m spatial resolution TM data. 

The additional spectral heterogeneity in the higher resolution data, produced 
higher variability or class overlap about class means, resulting in worsening 
cover class confusions. New classifiers or approaches should be investigated 
which take into account the increased spectral heterogeneity of higher resolution 
data. 

Analysis of transformed divergences for all possible 2-band and 3-band combinations 
indicated a general trend of including spectral bands from the four major 
spectral regions (visible, near-infrared, middle-infrared, and thermal -infrared) 
in order to maximize cover class discriminations. Assessment of principal 
components for TM and MSS data indicated approxiamtely two additional components 
of information in TM vs MSS. As evidenced by an increase in classification 
accuracy from 78. 7 % to 83.2%, the addition of the thermal -infrared band (TM6, 
10.4-12.5*m) to the other 6 channels in a classification provides for a significant 
improvement in cover class discrimination. 


COMPARATIVE TECHNIQUES USED TO EVALUATE THEMATIC MAPPER DATA 
POR LAND COVER CLASSIFICATION IN LOGAN COUNTY, WEST VIRGINIA 


H. W. B lodge t 

NASA/Goddard Space Flight Center 

J. 0. Brumfield 
Marshall University 

R. G. Witt 

NASA/Goddard Space Flight Center 
R. F. Marcell 

Computer Sciences Corporation 
Goddard Space Flight Center 


The Eastern Regional Remote Sensing Applications Center (ERRSAC) and Marshall 
University (West Virginia) are conducting a joint evaluation of Thematic Mapper 
(TM) data for mapping abandoned lands and stages of revegetation in conjunc- 
tion with the West Virginia Department of Natural Resources (WVDNR). The 
study area for the project is Logan County, West Virginia, a predominantly 
rural county that contains a variety of land cover types, including many active 
and abandoned surface mines. The objectives of the research are to identify 
and map active, abandoned, and reclaimed/revegetated mine areas, and to test 
a variety of digital data processing techniques. The following data extraction 
and data reduction procedures are being compared: c&nonical analysis, principal 

components analysis, and band selection. 

A subset of the 4 November 1982 data set (TM seven bands) corresponding to 
the central area of Logan County was created. This subarea was then subjected 
to the above tiansformatlons, before being input to an unsupervised clustering 
algorithm which generated classifications for ecch technique. These classi- 
fications were labelled with the aid of USGS 7.5' topographic quadrangles. 

The WVDNR is in the process of photointerpreting three of the quads (Holden, 
Henlawi'on, Clothier) to provide ground truth information for comparison. 

It is expected that preliminary reeults for the Holden quadrangle will 
indicate that the canonically transformed data set yields the highest accu- 
racies for surface mine categories when compared with the photointerpreted 
and digitized ground truth data from WVDNR. 


RELATIVE ACCURACY ASSESSMENT AND SPECTRAL 
COMPARISON OF LAND SAT -4 MSS AND TM 

DATA FOR LAND COVER INVENTORY 

E. M. Middleton 

NASA/Goddard Space Flight Canter 
Creenbelt, MD 20771 

Y. C. Lu 

Computer Sclencee Corporation 
Silver Spring, MD 20910 

R. G. Witt 

NASA/Goddard Space Flight Center 
Creenbelt, MD 20771 

R. Sekhon 

Computer Sciences Corporation 
Silver Spring, MD 20910 

A study was undertaken to compare digital data for the Washington, DC 
scene from: a) the Landsat-4 Multlspectral Scanner (MSS) acquired 
11/02/82 (ID: 40109-15140); b) the Landsat-4 Thematic Mapper (TM) 

(ID: 40109-15140), simultaneously acquired; and c) the Landsat-3 
MSS from 7/11/61 (ID: 22362-1^)4) . Level I mapping accuracy was 
determined by a per pixel comparison of each remotely sensed data set 
with digitized ground truth for four USGS 7.5 minute topographic 
quadrangle maps. The overall mapping accuracy was computed as the 
sum of all correctly classified pixels (for all cover types and 4 
quads) divided by the total number of pixels for the 4 quads. Pre- 
liminary results Indicate comparable overall mapping accuracies from L-3 
and L-4 MSSs. An apparent Improvement In mapping accuracy from L-4 
TM is being quantified. In addition, the statistical variation In the 
spectral data from the three sensors associated vith each cover type 
were also evaluated. 
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Research Plans and Panel Discussion 




Thematic Mapper Research Plans 
Mark Settle 

Office of Spac^ Science and Applications 
NASA Headquarters 
Washington, DC 


Thematic Mapper (TM) data collected during the first two and one-half years 
of Landsat 4 operations will be primarily used by NASA for purposes of 
research and development. NASA's RAD objectives during this period are: 

(1) to evaluate the quality of TM data from an engineering perspective; (2) 
to verify the results of airborne TM simulator studies conducted In the 
past; (3) to evaluate the utility of TM and MSS (Multlspectral Scanner) 
data for conventional remote sensing applications; and (4) to explore new 
avenues of scientific research that can be addressed for the first time 
with the TM. The quality of TM data Is currently being evaluated by 
thirty-two scientific Investigators who were selected by NASA to 
participate In the Landsat 4 Image Data Quality Analysis (LIOQA) program. 
This program was Initiated In 1982 and will be largely completed In 1984. 

At the same time, Investigators at various NASA field centers are 
attempting to verify the results of earlier experiments that were performed 
with airborne sensors that simulated the measurement capabilities of the 
TM. Airborne TM simulators developed by NASA have been used to conduct 
experimental multlspectral surveys of local test site areas from elevations 
of 8 to 15 kilometers. Orbital TM Imagery Is row being analyzed to verify 
that similar results can be achieved over broad areas based upon 
multlspectral measurements of surface radiance obtained at an altitude of 
705 kilometers. 

The third major element of NASA's TM research program Involves comparative 
studies of TM and MSS Imagery to evaluate the relative utility of each 
sensor system for conventional applications such as crop monitoring, 
lineament mapping, and regional watershed analysis. The measurement 
capabilities of the TM are considerably more sophisticated than those of 
the MSS. Comparative studies of both data sets will be conducted to 
determine how the improved measurement capabilities of the TM enhance the 
quality and quantity of physical Information that can be derived from 
orbital multlspectral imagery. This research will attempt to relate 
Incremental Improvements in the Information content of TM Imagery to 
specific improvements In resolution and sensitivity that were incorporated 
In the design of the TM sensor. The results of this research will assist 
NASA and others In evaluating the relative utility of further Improvements 
In the resolution and sensitivity of Imaging sensor systems. 

Finally, the successful launch of Landsat 4 represents the beginning of a 
new era In land remote sensing research. Analysis of MSS data collected by 
Landsats 1, 2, and 3 during the 1970 ; s revealed that unique types of useful 
Information concerning the surface conditions on our planet could be 
derived from orbital multlspectral surveys. The decade of the seventies 
was primarily a period of technique development that resulted In a wide 
variety of new methods for analyzing and Interpreting digital Earth 
Imagery. With the advent of the TM, we now have a unique opportunity to 
build upon the results of this earlier research and apply multlspectral 
Imaging methods to the study of a wide variety of terrestrial phenomena. 
During the decade of the 1980's these methods have matured to the point 
where they can be meaningfully employed to study surface conditions and 
processes on a global scale. NASA currently plans to enlarge the scope of 
ongoing TM Investigations and explore Innovative scientific applications of 
TM data during 1984 and 1985. 


